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CHAPTER  1 

TRAILING  EDGE  NOISE 
AT  LOW  MACH  NUMBERS 


SUMMARY 

A  review  is  made  of  the  diffraction  theory  of  the  trailing  edge  noise  generated  by  a 
flat-plate  airfoil  of  zero-thickness  and  non-compact  chord,  according  to  which  the  sound 
is  attributed  to  the  scattering  of  a  ‘frozen’  pattern  of  turbulence  wall  pressure  swept  over 
the  edge  in  the  mean  flow.  Extension  is  made  to  determine  the  sound  produced  by  very 
low  Mach  number  flow  over  the  edge  of  an  airfoil  of  finite  thickness.  In  applications  it  is 
desirable  to  represent  the  noise  in  terms  of  a  surface  integral  over  the  airfoil  involving  a 
Green’s  function  and  a  metric  of  the  edge  flow  that  can  be  calculated  locally  using  the 
equations  of  motion  of  an  incompressible  fluid.  It  is  argued  that  the  appropriate  metric 
for  a  rigid  airfoil  is  the  incompressible  ‘upwash’  velocity  (determined  by  the  Biot-Savart 
induction  formula  applied  to  the  boundary  layer  vorticity  outside  the  viscous  sublayer), 
and  not  the  surface  pressure.  Formulae  for  calculating  the  noise  are  given  when  the  airfoil 
thickness  is  acoustically  compact,  and  for  both  three  and  two-dimensional  edge  flows. 

The  theory  is  illustrated  by  a  detailed  discussion  of  a  two-dimensional  vortex  flow  over 
an  airfoil  with  a  rounded  trailing  edge.  The  problem  is  simple  enough  to  be  treated 
analytically,  yet  is  also  suitable  for  validating  computational  edge  noise  schemes. 
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1.  INTRODUCTION 


The  ‘self-noise’  generated  by  an  airfoil  in  a  nominally  steady,  high  Reynolds  number  flow 
is  attributed  to  the  instability  of  the  airfoil  boundary  layers  and  their  interactions  with  the 
trailing  edge  [1-3].  The  edge  is  usually  a  source  of  high  frequency  sound  associated  with 
smaller  scale  components  of  the  boundary  layer  turbulence.  Low  frequency  contributions 
from  a  trailing  edge,  that  may  in  practice  be  related  to  large  scale  vortical  structures  shed 
from  an  upstream  appendage,  are  small  because  the  upwash  velocity  they  produce  in  the 
neighborhood  of  the  edge  tends  to  be  cancelled  by  that  produced  by  vorticity  shed  from 
the  edge  [4,  5].  If  the  surface  S  of  the  airfoil  is  rigid,  and  is  at  rest  in  a  mean  stream,  the 
far  field  acoustic  pressure  p'(x,  t)  =  p{'x.,t)  -  Po  at  position  x  and  time  t  {p{x.,t)  being  the 
pressure  and  Po  its  mean  value  in  the  acoustic  far  field)  is  given  formally  by  Curie’s  [6,  7] 
formula 


dxidxi  L  47r|x  -  y|  dxi  i 


47r|x-y| 


dxi 


dSj(y) 
47r|x  -  y|' 


x|  — >  oo. 


(1) 


where  the  first  integral  is  over  the  volume  V  occupied  by  the  fluid,  the  surface  element  dSj 
is  directed  into  V, 

Pij  =  {p-Po)Sij-aij,  (2) 


and  Gij  is  the  viscous  stress  tensor.  The  square  brackets  [  ]  in  (1)  denote  evaluation  of  the 
contents  at  the  retarded  time  t  —  |x  —  y|/co,  where  Cq  is  the  speed  of  sound.  The  direct 
sound  generated  by  the  turbulence  quadrupoles  [8,  9]  is  represented  by  the  first  integral  in 
(1).  At  low  Mach  numbers  (e.g.  underwater)  Tij  ^  poViVj  per  unit  volume  (po  and  v  being 
respectively  the  mean  fluid  density  and  the  velocity),  and  the  quadrupole  noise  is  usually 
negligible  compared  with  that  from  the  edge.  The  latter  is  contained  in  the  surface  integral 
in  (1),  and  for  an  acoustically  compact  surface  (for  example,  an  airfoil  whose  chord  is  much 
smaller  than  the  typical  acoustic  wavelength),  the  edge  noise  is  equivalent  to  that  generated 
by  a  distribution  of  dipoles  on  S,  whose  strength  per  unit  area  is  the  unsteady  surface 
pressure  [8,  9].  In  that  case  the  ratio  of  the  acoustic  powers  generated  by  the  edge  and  the 
volume  quadrupoles  ~  0(1/M^),  where  M  ~  v/co  <C  1  is  the  characteristic  Mach  number 
[6].  At  high  frequencies,  when  it  is  not  permissible  to  assume  the  chord  to  be  compact,  the 
relative  efficiency  of  the  edge  noise  is  increased  to  0(1/M^)  [10]. 

For  the  compact  chord  airfoil  the  dipole  strengths  can  be  determined  to  a  sufficient 
approximation  for  use  in  (1)  by  a  preliminary  calculation  of  the  flow  near  the  edge  in  which 
the  influence  of  compressibility  is  ignored.  However,  although  (1)  is  exact  for  a  non-compact 
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airfoil,  it  is  not  permissible  to  neglect  compressibility  when  specifying  the  dipole  strengths 
in  the  surface  integral.  Therefore,  edge  noise  predictions  have  traditionally  been  made  by 
one  of  two  alternative  procedures.  In  the  first,  it  is  assumed  that  the  trailing  edge  is  well 
approximated  by  a  semi-infinite,  rigid  plate.  A  calculation  is  then  performed  in  which  the 
free  field  hydrodynamic  pressure  pi,  say,  generated  by  the  boundary  layer  turbulence  is 
diffracted  by  the  edge.  If  pi  were  known  exactly  this  procedure  would  yield  an  accurate 
prediction  of  the  edge  generated  sound.  Hitherto  it  has  not  been  possible  to  prescribe 
with  sufficient  accuracy  the  influence  of  the  edge  on  the  hydrodynamic  pressure,  and  it  has 
usually  been  assumed  that  the  turbulence  is  swept  past  the  edge  by  the  mean  flow  as  a 
frozen  distribution  of  vorticity  [1,  11  -  15].  In  the  second  method,  the  noise  is  calculated  by 
the  acoustic  analogy  theory  of  Lighthill  [8,  9]  using  a  compact  Green’s  function  tailored  to 
the  trailing  edge  geometry  [16,  17]  (as  opposed  to  the  free 'space  Green’s  function  used  in 
Curie’s  equation  (1)).  The  sound  can  then  be  expressed  directly  in  terms  of  the  vorticity  in 
the  edge  flow  (calculated  as  if  the  flow  is  incompressible);  this  approach  also  enables  the 
frozen  approximation  of  the  first  method  to  be  extended  to  more  complicated  trailing  edge 
geometries  [17]. 

Modern  computational  procedures  will  soon  yield  accurate  predictions  of  the  high 
Reynolds  number  hydrodynamic  motion  near  a  trailing  edge,  and  it  is  appropriate  to 
re-evaluate  the  kind  of  numerical  data  that  will  be  required  to  make  confident  predictions 
of  the  far  field  sound.  We  have  suggested  above  that  an  incompressible  approximation 
to  the  surface  pressure  is  not  sufficient  when  the  surface  is  not  compact.  This  is  because 
a  non-compact  body  extends  into  the  acoustic  far  field,  and  any  representation  of  the 
sound  as  a  Curie  (or  Kirchhoff)  type  of  surface  integral  must  ensure  that  the  boundary 
conditions  on  the  surface  continue  to  be  satisfied  in  the  far  field.  This  can  be  done  either 
by  using  a  Green’s  function  tailored  to  the  airfoil  characteristics  or  (as  in  (1))  by  using  the 
free  space  Green’s  function,  but  with  p[j  known  to  the  required  precision  in  the  acoustic 
domain.  In  the  latter  case  the  dipole  strength  must  be  prescribed  with  full  account  taken 
of  compressibility,  which  is  not  normally  possible  because  it  presupposes  a  knowledge  of  the 
acoustic  field  at  the  surface  [18].  In  fact,  the  dipole  strength  can  safely  be  estimated  from 
an  incompressible  edge  flow  model  only  when  the  airfoil  is  acoustically  compact.  In  general, 
the  surface  source  strength  turns  out  to  depend  on  the  upwash  produced  by  the  unsteady 
flow.  When  this  is  known  (in  an  incompressible  approximation)  the  edge  noise  can  be  found 
by  using  an  acoustic  Green’s  function  whose  normal  derivative  vanishes  on  the  rigid  surface 
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of  the  airfoil. 

In  this  chapter  the  thin-plate  diffraction  theory  of  trailing  edge  noise  is  first  reviewed  (§2), 
and  predictions  are  used  to  exhibit  explicitly  the  failure  of  approximations  based  on  (1). 
The  general  edge  noise  problem  at  low  Mach  numbers  is  then  formulated  in  terms  of  the 
theory  of  vortex  sound,  and  it  is  demonstrated  how  the  the  sound  can  be  determined  from 
an  incompressible  approximation  to  the  ‘upwash’  velocity  (§3).  The  theory  is  illustrated  by 
a  detailed  application  (§4)  to  determine  the  noise  produced  in  a  low  Mach  number  vortex 
flow  past  the  rounded  trailing  edge  of  a  thick  airfoil. 
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2.  THIN-PLATE  MODEL  OF  TRAILING  EDGE  NOISE 
2.1  Diffraction  theory  [11,  12] 

Consider  turbulent  trailing  edge  flow  over  the  upper  surface  X2  =  -1-0  of  the  rigid 
half-plane  xi  <0,  X2  =  0,  where  the  main  stream  outside  the  boundary  layers  has  low 
subsonic  speed  U  in  the  xi -direction  (Figure  1).  The  calculation  of  the  edge  noise  is 
formulated  as  a  scattering  or  diffraction  problem,  in  which  the  pressure  pi,  say,  that  would 
be  produced  by  the  same  turbulent  flow  if  the  surface  were  absent,  is  scattered  by  the  edge. 
The  scattered  pressure  p'  includes  both  acoustic  and  hydrodynamic  components,  the  latter 
accounting  for  the  modification  of  the  near  field  pressure  by  the  surface.  The  condition 
that  the  normal  velocity  vanishes  on  the  half-plane  is  taken  in  the  high  Reynolds  number 
approximation 

dpijdx2  +  dp'/dx2  =  0,  on  S,  (3) 

where  S  denotes  the  ‘upper’  and  ‘lower’  surfaces  {xi  <0,  X2  =  ±0).  In  turbulence-free 
regions,  and  when  the  mean  stream  Mach  number  M  =  U/co  is  very  small,  pressure 
fluctuations  p(x,6(;)  of  frequency  w  (with  suppressed  time  factor  satisfy  the  Helmholtz 
equation 

(V"  +  «>  =  0,  (4) 

where  Kq  —  uj/co  is  the  acoustic  wavenumber.  The  presence  of  the  boundary  layer  and 
the  turbulence  are  ignored  except  insofar  as  they  are  responsible  for  the  pressure  pi;  in 
particular  p'(x,a;)  is  assumed  to  satisfy  (3)  everywhere,  and  the  scattering  of  sound  by  the 
shear  flow  is  neglected.  This  approximation  is  not  valid  at  very  high  frequencies  when  the 
acoustic  wavelength  is  comparable  to  the  thickness  of  the  boundary  layer. 

The  pressure  pi(x,  w)  must  be  an  outgoing  solution  of  Helmholtz’s  equation  in  X2  <  0,  in 
the  region  ‘below’  the  boundary  layer  sources;  on  X2  =  0,  p\  is  equal  to  half  the  boundary 
layer  blocked  pressure  Ps  that  the  same  turbulence  would  exert  on  an  infinite  plane  wall  at 
X2  =  -1-0,  and  we  can  write 

1 

PiM  =  7,  I,  <  0,  k  =((:,,  0,  fc,),  (5) 

Z  J—oo 

where  j{k)  =  —  k"^  is  either  real  with  sign  sgn(ci;)  or  positive  imaginary.  The  problem 

of  calculating  p'  now  reduces  to  the  determination  of  the  scattered  pressure  (from  (3)  and 
(4))  produced  by  the  interaction  of  each  Fourier  component  ofpi  with 

S. 
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The  calculation  can  be  performed  by  the  Wiener-Hopf  procedure  [19],  which  supplies  the 
following  representation  of  the  total  perturbation  pressure  (which  is  finite  at  the  edge  of  the 
half-plane) 

p{x,u)  =  ]-  r 

2  7—00 

_  sgn^  p,(k,a;)y^(/c2-fc3)i/2^ - (6) 

47ri  y-oo  ^{Kl-kDy^  +  KiiKi  -kx+  iO) 

where  K  =  \JkI-\-  kj.  The  first  integral  represents  the  direct  pressure  pi  (generated  by  the 
boundary  layer  ‘in  the  absence  of  S’),  and  is  strictly  valid  only  outside  the  boundary  layer 
in  a;2  <  0  or  X2  >  S,  where  S  is  the  boundary  layer  thickness. 

The  integration  with  respect  to  Ki  can  be  performed  explicitly  when  xq  —>■  ±0.  It  is  zero 
in  the  wake  (xi  >  0),  where  the  scattered  pressure  vanishes.  For  xi  <  0  the  integration 
contour  is  displaced  to  — ioo  in  the  iCi-plane,  capturing  the  residue  contribution  from  the 
pole  at  Ki  —  ki  —  iO  and  an  integral  along  the  branch  cut  of  y  (Kq  —  kiyi"^  -I-  Ki,  which 
extends  from  — (k^  —  just  below  the  real  axis,  to  —ioo.  The  branch-cut  integral 

can  be  expressed  in  terms  of  the  error  function  erf  (a;)  =  (2/v^)  fo  [20].  The  total 

surface  pressure  is  then  found  to  be 

p{xi,±0,xz,uj)  =  ^  Jj  Ps(k,w)  l±erf(e"'?|a;i|5^(K2  -  A:|)i/2  +  A:i)]  e*’"'’"  d^k,  Xi  <  0. 

(7) 

The  argument  of  the  error  function  has  positive  real  part  for  all  real  values  of  ki,  so  that  the 
error  function  1  as  a;i  — >•  — oo.  Thus,  p  — >■  0  on  the  lower  surface  (a;2  =  —0)  far  upstream 
of  the  edge,  whereas  p  — >■  on  the  surface  X2  =  +0  exposed  to  the  turbulent  stream.  This 

occurs  at  distances  upstream  of  the  edge  exceeding  the  characteristic  eddy  dimension.  If 
the  impinging  boundary  layer  turbulence  is  assumed  to  be  frozen  during  convection  over  the 
edge,  measurements  of  the  blocked  pressure  Ps  several  boundary  layer  thicknesses  upstream 
of  the  edge  can  be  used  in  the  second  integral  of  (6)  to  predict  the  edge  noise. 

At  large  distances  from  the  edge  the  integrations  with  respect  to  K\  and  kz  in  the  second 
integral  of  (6)  may  be  performed  by  the  method  of  stationary  phase  [19,  21].  This  yields 
the  edge-scattered  acoustic  pressure  in  the  Chase-Chandiramani  [11,  12]  form 


P'(x,w) 


Ko  sin  2  ^sin(0/2)e*''‘’l' 

^/2\x\ 


y/KoSinip  -I-  ki 
(Koa:i/|x|  -  ki) 


Ps{ki,  Ko cosip,uj)dki, 
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where  the  angles  0,  defining  the  radiation  direction  x/|x|  are  indicated  in  Figure  1. 

The  integrand  in  this  formula  is  singular  at  k\  =  Koa;i/|x|,  where  the  stationary  phase 
approximation  breaks  down.  However,  this  occurs  when  ki  lies  in  the  acoustic  domain, 
i.e.,  for  a  blocked  pressure  Fourier  component  Ps{ki,KoCOsip,uj)  that  actually  represents 
a  sound  wave  generated  by  the  boundary  layer  quadrupoles  in  the  absence  of  the  edge. 
Such  contributions  can  be  neglected  at  small  Mach  numbers.  When  the  blocked  surface 
pressure  is  regarded  as  entirely  hydrodynamic,  the  remaining  integral  in  (8)  is  dominated  by 
wavenumbers  ki  u/U  ^  Ko>  and  (8)  reduces  to 


p'{x,u}) 


-«o  sin2 '0sin(^/2)e*''“l^l  r°°  ps{ki,0,uj)dki 

7—0 


V2\: 


^/ki  iO 


X  — )•  oo. 


(9) 


This  representation  is  suitable  for  expressing  of  the  edge  noise  in  terms  of  the  hydrodynamic 
(i.e.,  incompressible)  component  of  the  blocked  pressure  measured  upstream  of  the  edge. 


This  is  usually  done  by  referring  0,0;)  to  the  blocked  surface  pressure  wavenumber- 
frequency  spectrum  P(k,  cu)  [22].  It  is  assumed  that  a  finite  section  —^L  <  xz  <  \L,  say,  of 
the  trailing  edge  is  wetted  by  the  turbulent  flow,  where  L  is  much  larger  than  the  boundary 
layer  thickness  6.  Then,  for  statistically  stationary  turbulence 


(j),{kuO,u)p*{k[,0,oj'))  «  —6{u  -  oj')5{ki  -  k[)V{h,0,u),  L->5,  (10) 

where  the  angle  brackets  (  )  represent  an  ensemble  average,  and  the  asterisk  denotes 
complex  conjugate.  The  acoustic  pressure  frequency  spectrum  $(x,a;)  of  the  edge  noise 
(defined  such  that  (p'^(x,t))  =  $(x,a;)da;)  is  then  calculated  from  (9)  to  be  given  by 


$(x,a;) 


oiL  sin^(0/2)  sin'!/)  P{ki,0,u})  dk 


27rCo|xP 


£ 


\ki 


M  ^  1,  |x|  oo. 


(11) 


The  peak  acoustic  pressures  are  radiated  in  the  ‘forward’  directions  9  =  ±7r.  Numerical 
predictions  are  made  by  introducing  a  convenient  empirical  model  for  P(k,  w)  (see,  e.g.  [22, 
23]),  although  this  will  not  be  discussed  here. 


2.2  Application  of  Curie’s  equation 

To  derive  these  results  from  Curie’s  representation  (1),  the  first,  quadrupole,  integral  on 
the  right  is  discarded,  and  is  approximated  by  {p  —  Po)dij  at  high  Reynolds  numbers. 
Then  for  each  component  of  the  sound  of  frequency  u) 

Q  roo  rO  ei«o|x-y| 

y  =  (y„0.!,3)  (12) 


7 


Report  No.  AM  99-003 


Boston  University,  College  of  Engineering 


where 

\p{xi,X3)]  =  p{xu+0,xz,u)  -  p{Xu-0,X3,Uj)  (13) 


is  the  pressure  jump  across  the  half-plane. 
At  large  distances  from  the  wetted  edge 


p'{x,uj) 


/  giKo|x-y|  \ 

Ulx-ylj^ 

4,r|xP  ^ 

:  0  for  yi  >  0) 

sin  i/)  sin 

roo  roo 

/  dyz  /  \p{yi,y3)]e 

J —  00  J —00 

47r|x| 

—iriKo  sin  'ijjsind, 


(K;oa:i/|x|,Korc3/|x)|)]e’ 


tKo  X 


X  — >•  oo, 


(14) 


where 

1  rOO  rOO 

[p(^i,^3)]  =  7;^  /  dys  (15) 

(27rj  J —00  2—00 

The  first  line  of  (14)  implies  that  the  principal  contribution  to  the  integral  is  from  those 
components  of  [p(yi,y3)]  with  length  scales  ~  0(1/aco)-  In  other  words,  a  correct  evaluation 
of  the  integral  requires  the  retention  of  phase  information  in  [p]  characterizing  fluctuations 
in  the  surface  pressure  over  distances  of  the  order  of  the  acoustic  wavelength,  which  typically 
exceeds  the  scale  of  the  hydrodynamic  motions  by  a  factor  1/M  ^1.  In  this  integral 
phase  interference  with  the  exponential  factor  is  responsible  for  correcting  the 

directivity  of  the  sound  from  that  of  a  free  field  dipole  ~  sin  tpsind  =  cos  ©  (0  being  the 
angle  between  x/|x|  and  the  a:2-axis)  orientated  normal  to  the  airfoil,  to  sin^  ipsm{6/2), 
which  conforms  to  the  rigid  body  surface  condition  {dp' /dxn  —  0)  in  the  far  field. 

If  an  attempt  is  made  to  approximate  the  surface  pressure  jump  [p{yi,yz)]  in  Curie’s 
formula  (12)  by  its  value  for  incompressible  flow,  it  would  be  equivalent  to  setting 


[p(«;oa;i/|xl,«:oa;3/|x)|)]  «  [p(0,0)]  = 


roo  rO 


-  (2,rr 


in  (14),  where  F{uj)  is  the  net  normal  force  exerted  on  the  fluid  by  the  half-plane.  This 
force  may  be  determined  exactly  from  the  solution  (6)  to  be  given  by 


F(.>) = -2«  r 

J—00  Ko{ki  —  iO) 


dki- 


(16) 
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It  is  unbounded  in  the  incompressible  limit  in  which  «o  — >  0,  and  cannot  be  reliably 
computed  at  low  Mach  numbers  by  a  numerical  simulation  of  the  flow,  being  determined 
by  the  unsteady  surface  pressures  over  a  distance  from  the  edge  of  the  order  of  the  acoustic 
wavelength. 

A  correct  prediction  of  the  edge  noise  from  Curie’s  formula  is  possible  only  when  the 
surface  pressures  are  known  within  the  acoustic  domain.  To  use  (14),  [p(/Co3^i/|x|,  Koa;3/|x)|)] 
must  be  evaluated  from  the  exact  formula  (6),  in  the  form 


K  X  ’ 

X  jj 

-1 


/: 


°°  V Ko  sin  Ip  +  ki 


2-\/2/Co7ri  sin2  ipcos{p/2)  (ko3;i/|x|  —  ki) 
whose  use  in  (14)  leads  directly  to  (8). 


Ps{ki,KoX3/\xluj)dki, 

(17) 


2.3  Kirchhoflf  integral  representation 

The  scattered  pressure  p'(x,a;)  of  §2.1,  satisfies  the  Helmholtz  equation  (4)  everywhere. 
By  introducing  a  Green’s  function  G{x,y,u),  which  is  any  solution  of 

(V"  +  kI)G  =  i(x  -  y)  (18) 

with  outgoing  wave  behavior,  p'{x,uj)  may  be  represented  by  the  following  Kirchholf  integral 
over  S  [24] 

p'{x,uj)  =  ^  (G{x,y,u)^{y,oj)  -  p'{y,uj)^{x,y,uj)^  dS{y),  (19) 

where  the  normal  derivatives  d/dyn  are  directed  into  the  fluid.  This  equation  determines  the 
scattered  sound  provided  p'  and  dp'/dyn  are  known  on  S.  However,  for  an  arbitrary  choice 
of  the  Green’s  function  G{x,y,u),  and  for  the  reasons  discussed  above  for  Curie’s  equation, 
acceptable  predictions  of  the  sound  are  possible  only  if  the  variations  of  p'  and  dp'/dyn  on 
S  are  specified  correctly  over  length  scales  comparable  to  the  acoustic  wavelengths. 

For  example,  one  might  attempt  to  express  the  radiation  entirely  in  terms  of 
the  surface  pressure  (already  determined,  say,  by  means  of  a  subsidiary  calculation 
valid  in  the  neighborhood  of  the  edge)  by  using  a  ‘pressure  release’  Green’s  function 
that  vanishes  on  S.  The  first  term  in  the  integrand  of  (19)  is  then  absent,  and 
p'{x,(jo)  =  —  §^p'{y,u))dG{x,y,Lo)/dyndS{y)  would  be  an  exact  representation  of  the  sound. 
But,  if  p'{y,Lo)  in  the  integrand  is  known  only  in  an  incompressible  approximation,  the 
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predicted  behavior  of  p'{'x.,u))  at  large  distances  from  the  edge  would  be  governed  by  the 
behavior  of  G{x,y,u)  as  |x|  -)•  oo,  and  the  far  field  scattered  pressure  would  therefore  be 
predicted  to  vanish  on  S,  whereas  for  a  rigid  surface  |p'(x,a;)|  actually  assumes  its  largest 
values  there!  In  order  for  this  latter  behavior  to  be  predicted,  sufficient  phase  information, 
characterizing  changes  in  p'(y,  u)  on  S  over  distances  of  the  order  of  the  acoustic  wavelength, 
must  be  included  to  ensure  that  a  correct  estimate  is  obtained  for  the  asymptotic  behavior 
of  the  integral  (19)  as  |x|  — )■  oo.  This  requirement  is  equivalent  to  the  correction  of  the  free 
field  dipole  directivity  (cos0  =  simjjsinO)  of  (14)  brought  about  by  the  use  of  the  exact 
formula  (17)  for  [p]  in  (12). 

The  need  for  such  detailed  phase  information  in  the  prescribed  boundary  values  of  p' 
and  dp' /dyn  can  be  avoided  by  using  a  Green’s  function  that  already  satisfies  the  relevant 
boundary  conditions  on  S.  It  can  then  be  expected  that  surface  values  calculated  from  an 
incompressible  model  of  the  flow  will  be  sufficient  to  determine  the  far  field  sound  at  low 
Mach  numbers.  For  the  rigid  half-plane  the  Green’s  function  should  have  vanishing  normal 
derivative  on  S  (reciprocity  actually  implies  that  dG/dxn  =  0,  dG/dyn  —  0  respectively 
on  Xi  <  0,  X2  =  0  and  yi  <0,  y2  =  0).  Then  the  second  term  in  the  integrand  of  (19)  is 
absent,  and  condition  (3)  gives 

p'(x,w)  =  |^(yi, 0,2/3, w)G(x,y;a))dS(y). 

Js  dyn 

where 

[G(x,  y,  w)]  =  G(x,  2/1,  -fO,  2/3,  w)  -  G(x,  yi,  -0,  ya,  <^)  (21) 

is  the  jump  in  the  value  of  G  across  the  half-plane. 

The  integrals  in  (20)  may  be  evaluated  by  taking  dpi/dy2  to  have  its  value  when 
compressibility  is  ignored,  provided  that  those  turbulence  eddies  responsible  for 
the  edge  noise  are  always  very  much  closer  than  an  acoustic  wavelength  from  the 
edge,  which  is  always  the  case  at  sufficiently  small  Mach  numbers.  It  then  becomes 
appropriate  to  expand  Green’s  function  in  terms  of  the  nondimensional  source  distance 
/to\/yi  +  yl  (~  \/yi  +  yl/acoustic  wavelength)  from  the  edge.  When  the  observation  point 
X  is  in  the  acoustic  far  field  we  find  [25] 

G(x,y,w)  =  Go(x,y,a;)+Gi(x,y,a;)  +  ...,  (22) 
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where,  for  |x  -  ysUl  oo  and  /CoV^/i  +  2/2  < 
Go(x,y,u;)  = 


_ ^  piKoIx-yaUI  (7,  Cx  VO)')  —  _ - _ (y)  giKoIx-ysbl 

47rlx-j/3i3|  ’  |x-y3i3P/2 


i3  is  a  unit  vector  parallel  to  the  a;3-axis  (the  edge),  and  the  function 

99*  (x)  =  -\/rsin(0/2)  =  ^j^sin^  "0  sin  (0/2) 


(23) 


(24) 


is  equivalent  to  the  velocity  potential  of  incompressible  flow  around  the  edge  (in  the 
anticlockwise  direction)  expressed  in  terms  of  polar  coordinates  (a:i,a:2)  =  ^(cos^,  sin0). 
The  component  Go  represents  the  radiation  from  a  point  source  at  y  when  scattering  by  the 
half-plane  is  neglected.  The  component  Gi  provides  the  first  correction  due  to  the  presence 
of  the  half-plane,  and  (since  [Go]  =  0)  gives  the  leading  approximation  to  the  edge  noise 
when  used  in  (20). 

By  introducing  the  representation  (5)  of  pi  in  terms  of  the  blocked  pressure,  we  find, 
using  (23)  in  (20), 


p'(x,  w) 


^2/Co  sin 2  '0sin(0/2)e*''“l’'l 


TT 


Loo  Loo  ^  ^  ^ 


iKo  sin2  'i/)sin(0/2)e®'^"l’'l 


TT 


j{k)ps{k,u})d^k  / 

-OO  J—oo  J—00 

(25) 


The  t/3-integral  equals  2Tr5{k3)',  the  j/i-integral  must  be  treated  as  the  Fourier  transform  of 
a  generalized  function  [26],  and  evaluated  by  integration  by  parts,  when  it  is  found  to  be 
equal  to 

—1  /  ni 

2\ki  \  \j  ki  -\-  iO 

When  compressible  effects  in  the  specification  of  the  blocked  pressure  are  neglected,  it  may 
be  assumed  that  Kq  C  ki  for  all  relevant  values  of  ki  in  the  wavenumber  integral.  Then 
7(A:)  — >■  z|A:i|  and  (25)  reduces  to  precisely  the  low  Mach  number  approximation  (9)  obtained 
previously  by  diffraction  theory. 
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3.  ARBITRARY  TRAILING  EDGE  GEOMETRY 


3.1  Vortex  sound  theory 

The  diffraction  theory  of  §2.1  is  conveniently  extended  to  an  airfoil  of  finite  thickness 
(Figure  2)  by  means  of  the  theory  of  vortex  sound  [25],  in  which  the  total  enthalpy  B, 
rather  than  the  pressure  or  density,  is  taken  as  the  fundamental  acoustic  variable.  When 
the  mean  flow  Mach  number  is  small  enough  that  the  convection  of  sound  may  be  ignored, 
and  when  the  mean  fluid  density  is  constant. 


\c2  df^ 


B  =  div(S7  A  v). 


(26) 


where  f2(x,  t)  =  curl  v  is  the  vorticity.  In  those  regions  where  the  unsteady  motion  is 
irrotational  (D  =  0)  yet  still,  perhaps,  predominantly  nonlinear,  it  can  be  described  by  a 
velocity  potential  say,  which  satisfies  B  =  —d<f)/dt.  In  the  acoustic  far  field  the 

small  amplitude  pressure  fluctuations  are  determined  in  terms  of  B  by  the  linear  relation 


t) 

Po 


(27) 


Equation  (26)  relates  fluctuations  in  B  to  the  vorticity  and  velocity.  The  radiation 
condition  requires  the  solution  to  have  outgoing  wave  behavior,  and  for  each  Fourier 
component  of  frequency  u  it  can  be  expressed  as  the  sum  of  a  Kirchhoff  integral 
representing  a  contribution  from  the  surface  S  of  the  airfoil  (as  in  (19))  plus  the  direct 
radiation  from  the  vortex  sources: 

B{x,uj)  =  ^  (^^{y,uj)G{x,y,uj)  -  R(y,a;)^(x,y,a;)^  dS(y) 

-  j  G{x,y,u})(div{n  /\v){y,u)))d^y,  (28) 

where  G{x,y,uj)  is  an  outgoing  solution  of  (18). 

Let  the  Green’s  function  satisfy  the  rigid  surface  condition  dG/dyn  =  0  on  S, 
and  use  the  divergence  theorem  to  transform  the  volume  integral  in  (28)  as  follows: 

/  Gdiv(f2  A  v)dV  =  -  h  A  v)  •  ndS(y)  -  f(fiAv)-  VGd^y, 
where  the  unit  normal  n  on  S  is  directed  into  the  fluid.  Then 

B(x,u;)  =  ^G(x,y,a;)  +  (f2  Av)  •  n^(y,a;)  dS(y)  +  j  |^(x,y,a;)  •  Av)(y,a;)dV- 

^  (29) 
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Very  close  to  S  viscous  effects  are  dominated  by  shear  stresses,  and  the  momentum  equation 
can  be  taken  in  Crocco’s  form 

dv/dt  -f-  VB  —  —fl  A  V  —  1/  curl  Q 

where  u  is  the  kinematic  viscosity.  But  v  =  0  on  S  when  the  airfoil  is  rigid  and  at  rest. 
Hence  (29)  becomes 

r  dG  f 

B(x,  w)  =  y  5^’  w)  •  A  vj  (y,  uj)(fy  -  G{x,  y,  c;;)curl  n  •  ndS(y), 

and  the  identity  Gcurl  fl  =  curl  (Gfl)  -I-  Si  A  VG  and  the  divergence  theorem  then  yield 

B(x,w)  =  y  —  (x,y,a;)- (SlAv)(y,tc;)dV-  2/ j^fl{y,u)  A -^{x,y,u)  ■  ndS{y).  (30) 

The  surface  integral  is  the  contribution  to  the  radiation  from  the  unsteady  skin  friction  on 
S,  and  is  usually  ignored  when  the  Reynolds  number  is  large. 


The  remaining  integrals  can  be  evaluated  at  low  Mach  numbers  by  expanding  the  Green’s 
function  in  the  form  (22),  provided  the  characteristic  acoustic  wavelength  is  much  larger 
than  the  airfoil  thickness.  The  principal  contribution  is  from  Gi(x,  y,  w)  which,  however, 
must  now  be  taken  in  the  form  [17] 


G.(x  y  c.)  =  -1 


(31) 


where  the  potential  (/?*(y)  =  ^*(2/1, 2/2)  of  (23)  is  replaced  by  'l'*(y)  =  $*(2/1,2/212/3))  which 
describes  incompressible  potential  flow  around  the  edge  of  the  airfoil,  such  that 


$*(2/1, 2/2, 2/3)  ->  ¥>*{yuy2), 


00. 


We  now  find,  in  the  acoustic  far  field. 


(Sl  A  v)  {y,uj)d^y 


-  1/ j^Sl(y,w)  A  •ndS(y)|,  |x| 00.  (32) 

The  radiated  sound  automatically  satisfies  the  rigid  surface  boundary  condition  on  the 
distant  parts  of  the  airfoil,  so  that  the  vorticity  fl  and  the  velocity  v  in  the  integrands  can 
be  approximated  by  their  values  calculated  for  incompressible  flow  near  the  edge  of  the 
airfoil.  The  term  in  the  brace  brackets  of  (32)  is  proportional  to  the  net  normal  force  F{u) 
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exerted  on  the  fluid  by  the  airfoil.  As  noted  in  §2.2,  this  force  increases  in  proportion  to  the 
square  root  of  the  acoustic  wavelength;  in  the  present  notation  it  is  given  by  [17] 

F(u})  ^  j  ■  (nAv)(y,a;)dV-  i^£fi(y,a;)  A  ■  ndS(y)|. 


3.2  Diffraction  theory 


The  evaluation  of  (32)  requires  the  vorticity  and  velocity  fields  to  be  known  in  the 
neighborhood  of  the  edge.  An  alternative  representation  of  the  sound,  which  is  analogous 
to  (9)  for  the  flat-plate  airfoil,  can  be  derived  by  consideration  of  the  diffraction  theory  of 
§2.1.  We  shall  not,  however,  assume  the  boundary  layer  turbulence  to  be  frozen  during 
convection  past  the  edge,  but  will  introduce  an  ‘incident’  disturbance  Bi  which  is  defined 
to  be  the  exact  solution  of  equation  (26)  in  the  absence  of  the  airfoil  when  the  vorticity  Cl 
and  velocity  v  on  the  right  of  (26)  have  their  exact  values. 


To  calculate  Bj  the  interior  of  the  airfoil  is  imagined  to  be  replaced  by  fluid  with  no 
acoustic  sources,  wherein  the  actual  motion  is  determined  by  the  source  distribution  f2  A  v 
outside  S.  Then,  for  each  component  of  frequency  w. 


1  r  d  /  \ 


(33) 


where  V  is  the  fluid  volume  outside  S.  When  compressibility  is  neglected  in  the  source 
region  near  the  edge 

d'^v. 


dVidvj 

and  therefore,  because  v  vanishes  on  S,  the  acoustic  pressure  pi  ^  p^Bi  in  the  far  field 
corresponds  to  the  quadrupole  field 


(34) 


Pi(x,w) 


— ( XiXj  1 


47rlx|  V  jx] 


-  2%)  |x|  ^  00, 


(35) 


which  is  negligible  compared  to  the  edge  generated  sound. 


The  effect  of  the  airfoil  is  calculated  by  setting 


B{x,  t)  =  B'{x,  t)  -b  Bi{x,  t),  (36) 

where  B'  satisfies  the  homogeneous  form  of  (26)  (no  sources).  B'  and  Bi  are  related  by  the 
no-slip  boundary  condition  on  S,  which  the  momentum  equation  gives  in  the  form 

VB' +  VBi  = —ucnvlCl  on  S.  (37) 
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Thus,  when  KirchhofF’s  formula  (19)  is  used  (with  p'  replaced  by  B'),  we  find 
B'(x,  w)  =  -  ^  G(x,  y,  w)  (1^  +  i/curl  fl  ■  nVy,  w)dS(y), 


(38) 


provided  dG/dyn  =  0  on  S. 

Taking  the  low  Mach  number  expansion  of  the  Green’s  function,  the  leading  order  term 
of  which  is  given  by  (31),  we  find 


//  \  sin2  ^sin(0/2)e*''°l’'l  //;,=  */  \dBi.  ,  5$*(y)  \  . 


|x|  ^  oo, 


(39) 


where  the  integrand  is  to  be  evaluated  using  incompressible  approximations  for  dBi/dy„ 
and  fl. 


This  prediction  of  the  far  field  sound  is  identical  to  that  given  by  (32).  Indeed,  because 
d^*/dyn  =  0  on  S,  the  following  relations  are  readily  seen  to  transform  the  first  term  of  the 
integral  of  (39)  into  the  corresponding  term  in  (32): 


div($*VBi  -  BiV^*)  =  -  BiVH*  =  «  -$*div(n  A  v). 


where  the  final  approximation  follows  from  the  incompressible  limit  of  equation  (26). 

By  setting  Kq  =  0  in  (33)  we  derive  a  local  incompressible  representation  of  Bi,  from 
which  it  is  readily  deduced  that 


VSi  -f  f2  A  V  =  curl 


I 


curl  (i7  A  v)d^y 
,  — =  -curl 
v  47r  X  —  y| 


47r|x  -  y(’ 


(40) 


where  use  has  been  made  of  the  curl  of  the  incompressible  momentum  equation: 
dv/dt  -1-  VB  =  —Vt  A  V  —  i/curl  fi.  In  the  viscous  sublayer,  close  to  the  surface  of  the  airfoil, 
the  motion  becomes  linear  and 


^  «  0. 

Outside  the  sublayer  viscous  diffusion  is  negligible,  and  uV'^Vt  may  be 
integrand.  Hence,  we  can  introduce  an  ‘upwash  velocity’  vj  by  means 
formula  [27] 


vi(x,  t)  =  curl 


r  ^{y,t)(fy 

Vj  47r|x  -  y|  ’ 


discarded  from  the 
of  the  Biot-Savart 

(41) 
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where  the  integration  is  confined  to  the  boundary  layer  vorticity  in  the  nonlinear  region  Yg 
outside  the  viscous  sublayer.  On  S  and  within  the  volume  of  the  airfoil,  dvi/dt  =  —  VBi,  in 
terms  of  which  the  far  field  sound  becomes 


p'{x,oj) 


PoW-v/^sina 

7r\/^|x| 


f^^*{y)vin{y,(^)dS{y),  |x|-)-oo,  (42) 


where  vin  =  vi  •  n.  Note  that  in  applications  to  problems  in  which,  for  the  purpose  of 
calculation,  the  whole  motion  is  regarded  as  inviscid,  the  definition  (41)  makes  it  clear  that 
in  calculating  vi  the  bound  vorticity  on  S  must  be  excluded  from  the  integral. 

Equation  (42)  generalizes  the  Chase-Chandiramani  formula  (9)  for  the  flat  plate  airfoil, 
to  which  it  reduces  when  S  is  a  semi-infinite  half-plane.  To  see  this  it  is  necessary  to  note 
that,  in  an  incompressible  approximation  of  the  flow  near  the  edge,  the  incident  pressure  pi 
of  (5)  is  the  solution  of 


^  -P„div(fi  A  V)  -  (i.^) 


when  the  presence  of  the  airfoil  is  ignored.  Hence  Bi  =  pi/po  +  where  v  is  the  fluid 
velocity,  which  vanishes  on  and  within  S  where,  accordingly, 

^  =  -VB,  =  — Vp,.  (43) 

ut  Pq 


When  vin  in  (41)  is  replaced  by  {l/ipoio)dpi/dyn  the  first  line  of  (25)  is  recovered,  leading 
to  our  previous  result  (9). 

For  a  time-stationary  random  flow  past  the  edge  we  can  introduce  a  frequency  correlation 
function  'JZvviy,y',^)  that  satisfies 


(^^In(y,‘^)^^^„(y^^^')>  =  T^vv{y,y',i^)S{uj  -u'). 

The  acoustic  pressure  frequency  spectrum  (defined  as  for  (11))  then  becomes 
pIu^  simp  2) 


(44) 


$(x,u)) 


7r®Colx|' 


£'7^„„(y,y^a;)$*(y)$*(y')c^S(y)c^S(y'), 


a;  >  0. 


(45) 
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3.3  Two-dimensional  source  distributions 


We  record  here  the  modifications  of  formulae  given  above  when  the  aeroacoustic  sources 
and  trailing  edge  can  be  regarded  as  two-dimensional,  with  no  dependence  on  the  spanwise 
coordinate  2:3.  At  very  small  Mach  numbers,  the  dominant  component  of  the  Green’s 
function  whose  normal  derivative  vanishes  on  S  is  obtained  by  integrating  (31)  over 
-00  <  yz  <  00.  When  x  lies  in  the  acoustic  far  field  the  integration  can  be  performed  by 
the  method  of  stationary  phase,  which  yields  [5] 


Gi(x,y,a;) 


y*(x)$*(y)  ^  -sin(g/2)$*(y) 


TT  X 


TT  X  2 


(46) 


where  now  x  =  (2;i,a;2),  y  =  (yi,  Vt)- 

Thus,  the  two-dimensional  analog  of  (32)  is  the  cylindrically  spreading  wave  field 


p'(x,a;) 


— PoSin(0/2)e*''°l^l 

7r|x|5 


a$*(y) 

dy 


A  v)(y,n;)dV 


-  z/j(fi(y,a;)  •nds(y)|,  [x] -)•  00, 

where  ds  >  0  is  the  line  element  on  the  airfoil  profile  in  the  a:irc2-plane.  This  has  the  simple 
time-domain  form 


-poSin(d/2) 

7r|x|2 


r  5$*(y) 

J  dy 


(f2Av)(y,t)dV- 


•  n  ds(y)  , 


which  decays  in  amplitude  like  l/|x|5  at  large  distances  from  the  edge,  and  where  the  term 
in  the  square  brackets  [  ]  is  evaluated  at  the  retarded  time  t  —  |x|/co.  For  inviscid  (or  very 
high  Reynolds  number)  flow  we  can  take 


p'(x,  t)  « 


-po  sin(g/2) 

7r|x|5 


fd^*{y) 

J  dy 


Av)  {y,t)d^y 


(47) 


Similarly,  the  diffraction  theory  formulae  (38)  and  (42)  have  the  corresponding 
representations 


p'{x,t) 


Po  sin  (0/2) 


7r|x 

_  Po  sin  (0/2) 


P-h*iy) 

2  JS 


dB, 


dyr, 


(y.<) 


ds 


TT  X 


2  Js 


dVin 


dt 


{y,t) 


ds,  |x|  — >■  00, 


(48) 


where  the  skin  friction  contribution  has  been  discarded. 
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4.  ROUNDED  TRAILING  EDGE 

Rounded  or  bevelled  trailing  edges,  of  the  type  depicted  in  Figure  3a,  are  frequently 
used  in  experimental  studies  of  trailing  edge  noise  at  low  Mach  numbers  [1,  27,  28].  The 
simplified  geometry  is  suitable  for  validating  numerical  methods  for  edge  noise  prediction. 
When  the  Mach  number  is  small  enough  for  the  edge  flow  to  be  regarded  as  incompressible, 
it  is  necessary  to  be  able  to  determine  numerically  the  velocity  and  vorticity  distributions 
near  the  edge,  or  equivalently,  the  ‘upwash’  velocity  vjn-  In  this  section  the  procedure  is 
illustrated  for  a  highly  simplified  two-dimensional  edge  flow  that  can  be  treated  analytically. 


4.1  Formulation 

A  line  (or  ‘point’)  vortex  of  circulation  F  is  parallel  to  the  edge  of  the  airfoil  and  convects 
in  the  mean  flow.  Except  very  close  to  the  edge,  the  airfoil  has  uniform  thickness  h  and 
its  ‘upper’  and  ‘lower’  surfaces  coincide  with  the  planes  X2  =  The  upper  surface  is 

rounded  over  the  interval  —£<  xi  <  0,  and  meets  the  lower  surface  at  A:  xi  =0,  X2  —  —\h 
(Figure  3a).  At  large  distances  from  the  airfoil  the  mean  flow  is  at  speed  U  in  the 
a;i-direction,  and  the  mean  circulation  around  the  airfoil  is  assumed  to  be  adjusted  to  make 
A  a  stagnation  point. 


Viscosity  is  ignored  in  the  body  of  the  fluid,  so  that  the  two-dimensional  motion  of  the 
vortex  can  be  determined  by  mapping  the  fluid  region  bounded  by  the  airfoil  in  the  complex 
plane  of  ^  =  a;i  -I-  ix^  onto  the  upper  half  of  the  C-plane.  To  do  this  we  introduce  the 
complex  function 


2  TT  \  2  , 


j5  =  constant  >  —  1. 


(49) 


The  transformation 

|  =  /(C,/?) 

maps  the  exterior  of  the  semi-infinite  airfoil  of  uniform  rectangular  cross-section  and 
thickness  h  shown  in  Figure  3b  onto  Im  C  >  0,  such  that  the  boundary  points  C,  A,  B,  D 
on  S  correspond  respectively  to  the  points  C',  A',  B',  D'  on  the  real  ^-axis. 


It  may  then  be  verified  that  the  composite  transformation 

|  =  ^^^[/(C,l)  +  a/(C,/5)],  «>0,  P>+1, 


(50) 
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where  a  is  a  second  constant,  maps  the  real  C-axis  onto  an  airfoil  profile  S  whose  upper  and 
lower  surfaces  coincide  with  X2  =  respectively  for  ReC  <  — /?  and  Re<^  >  1.  The  interval 
-/?  <  Re  C  <  -1  maps  into  an  upper  ‘rounded’  section  of  the  trailing  edge  (as  in  Figure  3a) 
which  terminates  at  C  =  -1  at  the  top  of  a  vertical  end-face  (a;i  =  0,  -\h  <X2<  -|/id-  A) 
of  thickness  A,  that  corresponds  to  the  interval  —  1  <  Re  ^  <  1.  The  constants  a  and  P 
are  determined  by  the  prescribed  values  of  the  ratios  ^|h  and  A//i.  For  the  calculations 
presented  in  this  chapter  we  take 


a  =  600,  p  =  86.9370  (51) 

for  which  ^|h  =  4  and  A/fi  =  0.0074.  The  corresponding  airfoil  profile  is  that  shown  in 
Figure  3a  (because  A  <C  /i  the  end-face  cannot  be  distinguished). 

When  the  edge  A  is  a  stagnation  point,  the  complex  velocity  potential  of  mean  flow  is 
readily  confirmed  to  be  given  by 

-Uh  ,  Aa 
7r(l-Fa)  [  ^  {1  +  py 

which  becomes  asymptotically  Wm  ^Uz  when  \z\  ':$>  h.  The  vortex  F  is  convected  by  the 
mean  flow  and  by  the  induced  velocity  field  of  ‘image’  vortices  in  the  airfoil.  If  A  is  also 
a  stagnation  point  of  the  unsteady  flow,  i.e.,  if  the  unsteady  Kutta  condition  is  applied 
there,  additional  vorticity  is  shed  from  A,  and  will  also  influence  the  motion  of  F.  We  first 
consider  the  motion  and  sound  generation  in  the  absence  of  vortex  shedding. 


(c  -  i)t 


(52) 


4.2  No  vortex  shedding 

Let  Zo{t)  denote  the  complex  position  of  F  at  time  t.  The  velocity  potential  of  the 
(incompressible)  flow  produced  by  F  is 

™r  =  ^{lii(C-Co)-ln(C-C)}. 

where  Co  =  C(^o)  is  the  image  of  Zq  in  the  (^-plane,  and  the  asterisk  denotes  complex 
conjugate.  The  complex  potential  of  the  velocity  component  of  the  motion  of  F  produced 
by  image  vorticity  is  obtained  by  subtracting  the  ‘self  potential’  (— ir/27r)  ln(2:  —  Zo)  from 
Wr-  The  equation  of  motion  of  F  is  accordingly  obtained  in  the  form 

dz*  _  —  iF 
dt  27r 


C"  C 

So  So 


2^  Co  -  Co* 


dWrr, 

,  dz 


(53) 
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where  the  suffix  ‘o’  implies  evaluation  at  z  =  Zo,  and  the  primes  indicate  differentiation 
with  respect  to  2:.  This  equation  must  be  integrated  numerically  to  determine  the  path  of 
the  vortex,  and  it  is  convenient  to  do  this  in  the  C-plane,  where  it  is  equivalent  to 

dCo  -ire* 

dt  27r 

The  solid  curve  in  Figure  4a  depicts  the  calculated  path  when  the  vortex  is  released  at  a 
large  distance  upstream  of  the  edge,  above  the  airfoil  at  a  stand-off  distance  d=  \h  from 
the  upper  surface.  The  vortex  initially  translates  along  a  path  parallel  to  the  upper  flat 
surface  of  the  airfoil  at  speed  U  +  u,  where 

r 


The  results  discussed  below  are  obtained  for  u  =  — O.IC/.  This  corresponds  roughly  to  a 
typical  large  fluctuation  velocity  close  to  the  wall  of  a  turbulent  boundary  layer  [30].  The 
trajectory  shown  in  Figure  4a  was  computed  (using  a  fourth  order  Runge-Kutta  scheme 
[20])  by  adjusting  the  value  of  Co  far  upstream  of  the  edge  to  make  the  initial  stand-off 
distance  d  (determined  in  terms  of  Co  by  (50))  equal  to  ^h.  The  vortex  location  is  indicated 
in  the  figure  at  different  nondimensional  times  Ut/h,  measured  from  the  instant  that  it 
crosses  Xi  =  0. 

Two  other  trajectories  are  also  shown  in  the  figure.  The  ‘frozen’  path  corresponds  to  the 
Chase-Chandiramani  approximation  [11,  12],  in  which  turbulent  structures  are  assumed  to 
translate  past  the  edge  at  a  uniform  mean  convection  velocity  parallel  to  the  airfoil.  The 
position  of  the  vortex  on  this  path  has  been  calculated  by  taking  the  uniform  convection 
velocity  toheU  +  u  (not  U),  i.e.,  its  ‘exact’  value  when  located  above  the  flat  section  of  the 
surface  of  the  airfoil.  The  path  labeled  ‘rdf  is  determined  by  neglecting  the  influence  of 
image  vortices  on  the  motion  of  F,  whose  equation  of  motion  is  now  (53)  with  the  omission 
of  the  first  of  the  two  terms  on  the  right  hand  side.  This  corresponds  to  the  approximation 
of  ‘rapid  distortion  theory’,  where  each  turbulent  element  is  assumed  to  convect  across  an 
inhomogeneous  region  at  the  local  mean  stream  velocity  [31,  32];  for  the  purpose  of  this 
illustration  it  has  again  been  assumed  that  the  mean  velocity  far  from  the  edge  is  17  -f  u.  In 
all  cases,  therefore,  the  vortex  passes  by  the  edge  at  approximately  the  same  speed. 

Sound  is  generated  by  the  vortex  principally  during  its  passage  past  the  edge,  and  can  be 
calculated  from  either  of  the  formulae  (47),  (48),  by  noting  that  the  transformation  (50) 


Co  c 

[K  Co  -  CJ 


+ 


dWm 

.  dC 


1C 


M2 


(54) 
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implies  that 


r(y)=Re<^-C. 


h 


1  + 


4a 


\  (1  +  /5)V]’ 

where  (  is  the  image  oi  z  =  yi  4-  iy2,  and  y  =  (j/i,  ^2)- 
Now,  for  a  line  vortex, 


(56) 


n  =  ri3d(y  -  x„(f))  and  V  =  ^(<), 


(57) 


where  Xo(t)  is  the  vortex  location  at  time  t  calculated  from  (53),  (54),  and  we  therefore  find 
from  (47) 


p'{x,t) 


-PoTsm{9/2)  fix' fl  +  4a/{l  + P) 


2\  5 


3 

TTZ 


1  -f  a 


Im 


dCo 


dt 


x|  -)•  00,  (58) 


where  (^o  is  evaluated  at  the  retarded  time  [t]  =  t  —  |x|/co. 


The  ‘exact’  curve  in  Figure  4b  is  a  nondimensional  representation  of  the  pressure 
signature  plotted  against  U[t]/h.  Also  shown  are  the  corresponding  predictions  of  rapid 
distortion  theory  and  the  frozen  approximation.  The  high  frequency  components  of  the 
sound  are  generated  by  scattering  at  the  sharp  edge  of  the  airfoil,  and  the  peak  radiated 
pressure  occurs  in  the  neighborhood  of  [t]  =  0  when  the  vortex  passes  close  to  the  edge.  In 
the  frozen  approximation  the  contributions  to  the  sound  at  higher  frequencies  are  reduced 
because  small  scale  disturbances  generated  by  the  vortex,  that  are  responsible  for  the  high 
frequency  sound,  decay  rapidly  as  the  stand-off  distance  of  the  vortex  path  from  the  edge 
A  increases. 


The  sound  can  also  be  calculated  from  the  diffraction  theory  formula  (48).  To  do  this, 
recall  that  Bi  =  —dcpi/dt,  where  in  two-dimensions  the  ‘free  space’  velocity  potential  4>i  of 
the  vortex  is  given  by 


(^i(x,t)  =  Re  |^^ln(^  -  ^o(<))| ,  z  =  Xi4-  1x2. 


It  now  follows,  using  the  representation  (56)  of  $*(y),  that 

,  |x|  00,  (59) 

where  the  integration  is  in  the  anticlockwise  direction  about  the  contour  S  of  the  airfoil  in 
the  z-plane,  and  the  square  brackets  denote  evaluation  at  the  retarded  time  t  —  |x|/co. 


p'(x,  t)^- 


PoT  sin(^/2) 


27r2 


1  -b  4a/{l  +  py 
1  +  a 


Re 


dzo  f  C  dz 
dt  fs  {z  —  Zo)'^ 
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3 

The  integrand  ~  0(1/|2:|2)  as  \z\  — >  oo.  The  integral  may  therefore  be  evaluated  by 
residues  by  shifting  the  integration  contour  to  a  large  circle  at  infinity,  thereby  capturing 
a  contribution  from  the  pole  at  z  =  Zo{t).  This  procedure  yields  (58).  Alternatively,  (59) 
can  be  used  to  investigate  the  contributions  to  the  diffraction  radiation  from  the  lower 
and  upper  surfaces  of  the  airfoil,  which  are  determined  by  the  respective  contributions  to 
the  integral  from  the  lower  surface  between  2:  =  — oo  —  ih/2  and  2:  =  —ih/2  and  from  the 
upper,  rounded  surface  between  z  =  —ih/2  and  —00  -1-  i/i/2.  These  separate  integrals  have 
been  evaluated  numerically,  and  the  corresponding  pressures  are  labelled  ‘lower  surface’  and 
‘upper’  in  Figure  5  for  the  conditions  of  Figure  4,  for  the  case  in  which  the  vortex  moves 
along  the  solid  (‘exact’)  path  of  Figure  4a.  The  net  acoustic  pressure  is  the  algebraic  sum 
of  these  separate  contributions. 

4.3  Influence  of  vortex  shedding 

When  the  Kutta  condition  is  imposed  at  the  edge  A  of  the  airfoil  during  the  passage  of  F, 
vorticity  is  shed  into  the  flow  and  swept  downstream.  The  acoustic  radiation  consists  of  the 
direct  sound  generated  by  F,  considered  in  §4.2,  and  the  sound  generated  by  the  wake.  If 
«  <C  17,  where  u  =  F/47rd  is  the  characteristic  induced  velocity  of  the  vortex,  the  trajectory 
and  the  acoustic  pressure  signature  of  F  are  to  a  good  approximation  the  same  as  when  the 
vortex  convects  at  the  local  mean  stream  velocity  (according  to  ‘rapid  distortion  theory’, 
see  Figure  4).  We  shall  therefore  adopt  this  approximation  to  examine  the  influence  of 
the  wake,  by  assuming  that  both  F  and  the  shed  vorticity  convect  at  the  local  free  stream 
velocity  along  undisturbed  streamlines  of  the  mean  flow.  This  is  equivalent  to  the  linearized 
approximation  of  unsteady  thin  airfoil  theory  [33,  34],  where  the  airfoil  is  modelled  as  a 
plate  of  infinitesimal  thickness  parallel  to  the  mean  flow  (as  in  §2),  and  all  perturbation 
quantities  are  proportional  to  the  amplitude  of  an  incident  ‘gust’.  In  that  case,  however, 
both  the  gust  and  the  wake  vorticity  convect  parallel  to  the  plate  at  the  same,  uniform 
mean  stream  velocity,  and  the  acoustic  pressure  generated  by  the  wake  turns  out  to  be 
equal  and  opposite  to  that  produced  by  the  gust,  so  that  there  is  no  net  radiation  from  the 
edge  [4]. 

The  shedding  does  not  become  significant  until  F  is  close  to  the  edge.  It  will  be  modelled 
by  releasing  from  the  edge  A  a  succession  of  elementary  line  vortices  of  circulation  F^  at 
discrete  times  t  =  4)  k  >1.  Let  Zk,  Cfc  denote  respectively  the  position  of  Fjt  in  the  z-plane 
and  its  image  in  the  upper  half  of  the  C-plane  at  time  t.  These  vortices  lie  on  the  stagnation 
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streamline  emanating  from  A  (see  Figure  6a),  along  which  they  convect  at  the  local  mean 
velocity;  their  images  Cit  are  located  on  the  line  Re  C  =  1  extending  from  the  real  axis  into 
the  upper  half-plane.  At  time  tj,  when  there  are  j  shed  vortices  in  the  flow,  the  complex 
potential  of  the  unsteady  component  of  the  flow  is 


w  =  ln(C  -  C.)  -  ln(C  -  C)}  +  i:  ln(C  -  G)  -  ln(C  -  «)}•  (60) 


27r 


k-i 


-i£k 

277 


The  Kutta  condition  is  satisfied  at  A  by  setting  dw/dC,  —  0  at  =  1.  This  determines  the 
strength  Tj  of  the  yth  vortex  in  terms  of  F  and  all  of  the  previously  shed  vortices  according 
to 


,  ,  ,  rimCo 


(61) 


|1  -  Col^  ’  jfe!  Im  Cfc 

To  apply  this  formula  it  is  necessary  to  specify  the  initial  position  (j-  =  1  -t-  ie  of  the  jth 
vortex  on  the  stagnation  streamline  Re  C  =  Ij  where  e  is  small  and  positive,  whose  precise 
value  does  not  materially  affect  the  results.  At  the  next  time  step  in  the  calculation  all 
of  the  shed  vortices  will  have  moved  along  this  streamline  by  distances  determined  by  the 
mean  velocity  potential  (52). 


The  results  of  such  a  calculation  are  illustrated  in  Figure  6a,  for  the  rounded  airfoil 
considered  previously.  The  initial  stand-off  distance  d  of  F  far  upstream  of  the  edge  is  equal 
to  |/i,  as  before,  and  the  characteristic  velocity  u  of  (55)  is  again  taken  to  be  —O.IU.  The 
trajectories  of  the  incident  and  shed  vortices  are  the  streamlines  of  the  mean  flow  shown 
in  Figure  6a.  Time  is  measured  from  the  instant  that  F  crosses  Xi  =  0,  and  the  figure 
shows  the  position  of  F  at  various  times,  and  also  the  corresponding  location  of  the  peak 
shed  vorticity.  This  peak  is  shed  into  the  flow  when  F  is  close  to  the  edge,  and  translates 
downstream  with  F  on  a  parallel  path  at  roughly  the  same  velocity.  Figure  6b  indicates 
how  the  overall  circulation  of  the  wake  vorticity  is  opposite  in  sign  to  F  and  has  a  final 
magnitude  equal  to  about  80%  of  F.  The  slope  of  the  wake  circulation  curve  is  always 
negative,  showing  that  the  sign  of  each  elementary  vortex  Vj  is  always  opposite  to  F. 

The  sound  generated  by  the  impinging  vortex  and  the  wake  at  the  retarded  time 
tj  =  t—  |x|/co,  just  after  the  release  of  the  jth  vortex  from  A,  can  be  calculated  from  either 
of  the  following  generalizations  of  (58)  and  (59) 


23 


Report  No.  AM  99-003 


Boston  University,  College  of  Engineering 


p'(x,  t) 


X  00. 


(62) 


The  pressure  signature  is  plotted  against  U[t]/h  in  Figure  7.  Also  plotted  are  the 
separate  contributions  from  F  and  from  the  wake,  which  are  of  comparable  magnitudes,  but 
opposite  in  sign.  Because  of  the  progressive  increase  in  the  total  wake  circulation,  when 
U[t]/h  exceeds  about  5  the  separate  acoustic  pressures  attributable  to  F  and  the  wake  are 
effectively  equal  an  opposite.  Thus,  the  net  radiation  is  significantly  different  from  zero 
only  for  retarded  times  U[t]/h  ~  0(1)  when  F  is  very  close  to  the  sharp  edge  A.  This  may 
be  contrasted  with  the  analogous  result  for  an  airfoil  of  infinitesimal  thickness  [4] ,  for  which 
the  predicted  radiations  from  F  and  from  the  wake  are  equal  an  opposite  for  all  times,  and 
linear  theory  accordingly  predicts  that  no  sound  is  produced  at  the  trailing  edge. 
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5.  CONCLUSION 

In  this  chaapter  the  Chase-Chandiramani  diffraction  theory  for  estimating  trailing  edge 
noise  from  a  flat  plate,  zero  thickness  airfoil  has  been  extended  to  low  Mach  number  flows 
past  a  non-compact  airfoil  of  finite  thickness  and  arbitrary  trailing  edge  geometry.  For  the 
flat-plate,  the  radiation  can  be  approximated  by  considering  the  diffraction  at  the  edge 
of  the  boundary  layer  blocked  pressure,  which  is  assumed  to  convect  in  a  frozen  pattern 
past  the  edge.  The  same  approximation  for  a  thick  airfoil  significantly  underpredicts  the 
high  frequency  components  of  the  sound.  In  this  case  the  problem  must  be  formulated 
in  terms  of  the  diffraction  of  the  the  boundary  layer  ‘upwash’  velocity.  Both  approaches 
are  equivalent  for  the  flat  plate  airfoil,  but  the  extension  permits  account  to  be  taken  of 
modifications  of  the  turbulence  during  convection  past  the  variable  geometry  edge. 

In  applications  it  is  desirable  to  be  able  to  make  accurate  predictions  of  low  Mach  number 
edge  noise  by  first  performing  numerical  simulations  of  the  edge  flow  based  on  the  equations 
for  an  incompressible  fluid,  and  then  inserting  the  results  into  a  suitable  surface  integral 
representation  of  the  radiated  sound.  When  the  airfoil  chord  is  not  acoustically  compact, 
it  is  not  possible  to  make  such  predictions  solely  from  a  knowledge  of  the  incompressible 
approximation  to  the  unsteady  surface  pressure.  This  is  because  the  airfoil  itself  extends 
into  the  acoustic  far  field,  and  the  prescribed  surface  pressure  must  contain  sufficient 
acoustic  information  to  ensure  that  the  calculated  radiation  satisfies  the  appropriate 
dynamical  boundary  conditions  on  the  airfoil.  In  order  to  satisfy  these  conditions  using 
incompressible  data,  the  surface  integral  should  involve  an  acoustic  Green’s  function  that  is 
specifically  tailored  to  the  boundary  conditions.  For  a  rigid  airfoil  the  normal  derivative  of 
the  Green’s  function  should  vanish  on  its  surface,  and  the  incompressible  data  required  to 
determine  the  far  field  sound  is  the  boundary  layer  ‘upwash’  velocity.  This  velocity  is  equal 
to  that  given  by  the  (free-field)  Biot-Savart  formula  applied  to  the  boundary  layer  vorticity 
lying  outside  the  viscous  sublayer. 
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Figure  1.  Coordinates  for  trailing  edge  noise. 
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Figure  2.  Low  Mach  number  turbulent  flow  near  the  trailing  edge  of  an  airfoil  of 
finite  thickness. 
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Figure  3.  (a)  Airfoil  of  thickness  h  with  a  rounded  trailing  edge  section  of  length  t 

(b)  Mapping  the  region  outside  a  rectangular  airfoil  onto  the  upper  half  of  the 
(^-plane. 
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Figure  4.  (a)  Trajectories  of  the  vortex  F  past  the  trailing  edge  when  u  —  — O.lf/,  i/h  =  4, 
d//i=  0.5.  (b)  Acoustic  pressure  signatures  p/{po(r/h)^  sin(0/2)^/i/|7rx|/8}. 
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Figure  5.  The  acoustic  pressure  p/{/)o(r//r)^  sin(^/2)y'^/i/|7rx|/8}  and  the  separate 
contributions  from  the  upper  and  lower  surfaces  of  the  airfoil  of  Figure 
4  when  u  =  — O.IC/,  ^|h  =  4,  d/h  =  0.5  and  when  the  vortex  moves  along 
the  solid  (‘exact’)  trajectory  of  Figure  4a. 
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Ut/h 


Figure  6.  (a)  Trajectories  of  the  vortex  F  and  shed  vorticity  in  the  approximation 
of  rapid  distortion  theory  for  u  =  — O.IC/,  t/h  =  4,  d/h  =  0.5 
(b)  Total  wake  circulation  as  a  function  of  time. 
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Figure  7.  The  acoustic  pressure  p/{po(r/h)^sin(^/2)^h/|7rx|/8}  and  the  separate  con¬ 
tributions  from  the  vortex  F  and  the  vortex  wake  in  the  rapid  distortion  ap¬ 
proximation  of  Figure  6,  when  u  =  —O.IU,  i/h  =  4,  d/h  =  0.5. 
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CHAPTER  2 

ATTACHED  AND  SEPARATED  EDGE  FLOWS 


SUMMARY 

The  diffraction  theory  of  Chapter  1  is  applied  to  investigate  the  noise  produced  by 
turbulent  flow  over  an  edge  whose  upper  surface  profile  (the  suction  side)  is  rounded.  The 
sound  is  expressed  in  terms  of  the  ‘upwash’  velocity  fluctuations  that  the  same  boundary 
layer  turbulence  would  generate  if  the  airfoil  were  absent.  An  approximate  method  is 
proposed  for  expressing  these  fluctuations  in  terms  of  local  properties  of  the  blocked  pressure 
generated  on  the  surface  exposed  to  the  turbulent  flow.  Predictions  are  made  of  the  edge 
noise  spectrum  for  both  fully  attached  flow  (that  remains  attached  right  up  to  the  trailing 
edge)  and  for  cases  where  separation  occurs  on  the  rounded  profile.  When  premature 
separation  occurs  the  amplitude  of  the  edge-generated  sound  decreases  exponentially  with 
increasing  frequency,  and  predicted  edge  noise  levels  are  significantly  smaller  than  estimates 
obtained  when  the  airfoil  is  modeled  by  a  rigid  half  plane.  For  attached  flow  turbulence  on 
the  suction  side  of  the  airfoil  always  passes  close  to  the  edge  and  interacts  strongly  with 
it,  but  contributions  from  the  interaction  of  the  same  turbulence  with  the  pressure  side  of 
the  airfoil  are  reduced  because  of  the  finite  airfoil  thickness.  In  this  case  sound  levels  fall 
short  of  those  for  a  rigid  half  plane  typically  by  about  5  to  10  dB,  the  precise  values  being 
dependent  on  frequency,  and  on  the  ratio  of  the  boundary  layer  thickness  to  the  mean  airfoil 
thickness. 
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1.  INTRODUCTION 


Consider  turbulent  flow  at  very  small  Mach  number  over  the  upper  surface  of  the  trailing 
edge  an  uncambered  airfoil  with  a  rounded  upper  surface  edge  profile,  as  illustrated  in 
Figure  1.  The  profile  is  the  same  for  all  spanwise  locations,  so  that  the  edge  is  parallel 
to  the  airfoil  span.  Airfoils  of  this  type  are  frequently  used  in  fiow  noise  tests  [1-3],  and 
the  simplified  edge  geometry  is  convenient  for  validating  numerical  codes  for  edge  noise 
prediction.  Let  the  main  stream  outside  the  airfoil  boundary  layers  have  low  subsonic  speed 
U  in  the  a;i-direction  of  the  rectangular  coordinate  axes  {xi,X2,X3).  The  coordinate  origin 
0  lies  on  the  midplane  of  the  airfoil  in  the  ‘vertical’  plane  of  the  trailing  edge  A  (see  Figure 
2a,  below).  The  upper,  rounded  edge  region  extends  over  the  interval  —£<  Xi  <  0,  and 
meets  the  lower  surface  at  A:  Xi  =0,  X2  =  —5/1.  Upstream  of  the  edge,  for  Xi  <  the 
airfoil  has  uniform  thickness  h  and  its  upper  and  lower  surfaces  coincide  with  the  planes 
X2  =  ±^h.  The  mean  circulation  around  the  airfoil  is  assumed  to  be  adjusted  to  make  A  a 
stagnation  point. 

The  theory  of  the  self  noise  generated  at  a  trailing  edge  of  finite  thickness  is  discussed 
in  Chapter  1.  When  the  airfoil  thickness  h  is  acoustically  compact,  it  was  shown  that  the 
acoustic  pressure  p(x,  t)  in  the  far  at  x  at  time  t  can  be  rendered  in  the  form 

/OO 

p{x,u)e~^^^  du,  (1) 

-00 


p(x,a;) 


PoKo  sin 2  ^sin(0/2)e*''°l- 


£  I  iw$*(y)ui„(y, uj)  +  i^D(y, u)  A  .  n  I  dS(y), 


TT 


y/27Ti\ 


X 


dy 


|x|  — >■  00, 


(2) 


where  the  integration  is  over  the  upper  and  lower  surfaces  S  of  the  airfoil.  The  terms  in  this 
formula  are  defined  as  follows 


po  =  mean  fiuid  density, 

Ko  =  iofco  =  acoustic  wavenumber, 

Co  =  speed  of  sound, 

V  —  kinematic  viscosity,  > 

f2  =  curl  v  =  vorticity,  v  =  velocity, 
sin^  =  r/|x|,  r  =  sjx\  +  x^, 


(3) 
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Thus,  {xi,X2)  =  r{cos9,sm9),  and  ip  is  the  angle  between  the  radiation  direction  x/|x|  and 
the  a;3-axis  (the  airfoil  span,  out  of  the  plane  of  the  paper  in  Figure  1). 

#*(x)  =  $*(a;i,a;2)  is  an  auxiliary  function  that  depends  on  the  shape  of  the  airfoil  trailing 
edge  profile.  It  satisfies  Laplace’s  equation  and  describes  an  incompressible,  irrotational 
flow  around  the  edge  of  the  airfoil  (in  planes  parallel  to  the  a;ia;2-plane),  and  is  normalized 
such  that 

^*{xi,X2) (p*{xi,X2)  = 'yrsm{9/2),  for  »  h,  (4) 

where  (p*{xi,X2)  is  the  corresponding  potential  of  irrotational  flow  (in  the  anticlockwise 
direction)  about  the  edge  of  the  half  plane  Xi  <0,  X2  —  0. 

The  unit  normal  on  S  directed  into  the  fluid  is  denoted  by  n,  and  vin  =  vj  •  n  is  the 
normal  component  of  the  ‘upwash  velocity’.  This  is  defined  by 


dvi 

-w  = 


curl  (f2  A  v)dPy 
V  47r|x  — y| 


L 


=  curl 


dt 


47r|x  -  y|’ 


(5) 


where  V  is  the  region  occupied  by  the  fluid.  The  vortical  boundary  layer  motion  is 
assumed  to  be  well  approximated  at  low  Mach  numbers  by  the  Navier-Stokes  equation  for 
incompressible  flow.  In  the  viscous  sublayer,  close  to  the  surface  of  the  airfoil,  the  motion 
becomes  linear  and 


dft 

at 


—  «  0. 


Outside  the  sublayer  viscous  diffusion  is  negligible,  and  uV^ft  may  be  discarded  from  the 
second  integrand  of  (5).  The  upwash  velocity  is  then  given  by  the  Biot-Savart  formula  [4] 


v,(x,t)  =  curl/ 

Jvs  47r|x-y| 


(6) 


where  the  integration  is  confined  to  the  nonlinear  region  of  the  boundary  layer  Vj  outside 
the  viscous  sublayer.  This  result  implies  that,  in  applications  to  flows  of  very  large  Reynolds 
numbers,  where  the  dominant  turbulent  motions  may  be  regarded  as  essentially  inviscid, 
possible  contributions  to  the  integral  (6)  from  potential  flow  bound  vorticity  on  S  must  be 
excluded. 


Equation  (6)  defines  the  velocity  vj  in  terms  of  the  boundary  layer  vorticity  when  the 
presence  of  the  airfoil  is  ignored,  i.e.,  when  the  volume  enclosed  by  S  is  imagined  to  contain 
fluid  moving  irrotationally  at  the  induced  velocity  of  the  boundary  layer  vorticity.  Of 
course,  prior  to  evaluating  the  upwash  velocity  from  either  of  the  representations  (5)  or  (6), 
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the  vorticity  and  flow  velocity  v(x,t)  must  first  be  determined 'with  full  account 

taken  of  the  presence  of  the  airfoil  including,  of  course,  the  no-slip  condition  on  S. 

The  first  term  in  the  brace  brackets  of  the  integrand  of  equation  (2)  represents  the 
contribution  to  the  edge  noise  from  the  normal  surface  stresses  on  S;  the  second  term 
determines  the  net  contribution  from  the  surface  shear  stress,  and  is  usually  neglected  at 
very  high  Reynolds  numbers. 

In  Chapter  1  the  two-dimensional  version  of  equation  (2)  was  used  to  examine  the  sound 
generated  during  the  mean  flow  convection  of  a  rectilinear  vortex  past  the  trailing  edge 
of  Figure  1,  including  the  influence  of  vortex  shedding  induced  from  the  edge  A.  In  this 
chapter  equation  (2)  will  be  applied  to  estimate  the  acoustic  pressure  frequency  spectrum  of 
the  high  frequency  self  noise  generated  by  turbulent  flow  over  the  rounded  edge  profile.  The 
theory  is  formulated  (in  §2)  to  enable  account  to  be  taken  of  changes  in  the  mean  turbulent 
boundary  layer  characteristics  with  distance  from  the  edge.  In  the  theory  of  edge  noise 
generated  by  flow  over  the  edge  of  a  thin  plate  airfoil  [5-7],  the  properties  of  the  boundary 
layer  turbulence  are  expressed  in  terms  of  the  blocked  wall  pressure  wavenumber-frequency 
spectrum,  measured  or  specified  empirically  just  upstream  of  the  edge.  In  the  present 
formulation  (§3)  it  is  assumed  that  the  relevant  frequencies  are  sufficiently  high  that  the 
wavenumber-frequency  spectrum  of  the  upwash  velocity  can  be  represented  locally  in 
terms  of  the  blocked  pressure  spectrum,  thereby  permitting  acoustic  predictions  to  be  made 
in  terms  of  measured  or  calculated  turbulence  parameters.  The  Mach  number  is  sufficiently 
small  that  the  blocked  pressure  can  be  calculated  from  the  incompressible  equations  of 
motion.  Numerical  results  are  presented  in  §3  for  fully  attached  and  separated  edge  flows, 
including  a  comparison  with  predictions  of  the  frozen  approximation  of  [5  -  7] ,  according  to 
which  the  boundary  layer  turbulence  is  convected  over  the  edge  along  a  path  parallel  to  the 
undisturbed  mean  stream. 
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2.  FORMAL  REPRESENTATION  OF  THE  EDGE  NOISE 


Introduce  a  curvilinear  coordinate  system  (s,sx,a;3),  where  s  is  measured  along  the 
curvilinear  stream  lines  of  the  potential  function  $*(a;i,a:2)  in  the  clockwise  direction,  sx 
lies  in  the  a;i2;2-plane,  normal  to  the  streamlines  and  directed  away  from  the  profile  S  of  the 
airfoil,  on  which  we  can  take  sx  =  0  (Figure  2a).  Then,  Sx  >  0  in  the  fluid,  and  on  S 


vin  =  lim 
— 0 


ds± 


{s,S±,X3,U)) 


lim 

sx->-~0 


dks, 


(7) 


where  $(s,  sx,a:3,a;),  ^{s,s±_,k3,u)  are  respectively  the  upwash  velocity  potential  and  its 
Fourier  transform  with  respect  to  X3.  This  potential  is  well  defined  throughout  the  whole  of 
space  where  the  vorticity  S7  =  0,  including  the  region  occupied  by  the  airfoil. 

When  the  final  representation  in  (7)  is  substituted  for  in  equation  (2),  with  dS  =  dsdy3, 
the  integration  with  respect  to  7/3  yields  the  delta  function  S{k3).  In  an  exact  integral 
representation  of  the  edge  generated  sound  (involving  the  specification  of  compressible  data 
on  S)  the  surface  integral  would  actually  yield  6{k3  —  Kocosip),  which  merely  implies  that 
only  those  Fourier  components  of  Vin{s,X3,u))  that  have  supersonic  spanwise  phase  velocity 
\u)/k3\  >  Co  are  responsible  for  the  radiation  of  sound  from  the  edge.  In  other  words,  the 
only  relevant  wavenumbers  are  confined  to  the  acoustic  domain  |fc3|  <  |/Co|,  and  correspond 
to  variations  in  the  spanwise  direction  occurring  over  distances  of  the  order  of  the  acoustic 
wavelength.  However,  the  dominant  hydrodynamic  motions  near  the  edge  have  length 
scales  of  the  order  of  the  boundary  layer  displacement  thickness  5*  I/kq  <  l/l^sl-  In  the 
approximation  of  equation  (2),  KoS*  is  so  small  that  we  can  take  ^3  =  0. 

When  k3  =  0  and  the  source  motion  is  regarded  as  incompressible,  the  velocity  potential 
^{s,s±,k3,u})  satisfies  Laplace’s  equation,  which  locally  reduces  to 

f  d‘^\  - 


The  corresponding  stream  function  ^{s,s±,iv),  say,  is  related  to  l'(s,  sxjOjO;)  by  the 
Cauchy-Riemann  equations  [4] 


9#  . 

g^  =  -(5.5x.O.^). 


(8) 


Hence,  substituting  from  (7)  into  equation  (2),  and  discarding  the  contribution  to  the 
radiation  from  the  surface  shear  stress,  we  find 
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p{x,uj) 


PqUJ  l2iKo  ._1 
Ixl 


—  sin2 f  ^*(y(5))(-^ — (5,  5_l,  0,a;)  )  ds 
,  ,  TT  7-00  \ds^ 

-|^\/^^sin5'0sin(0/2)e*''‘>l’'l  f  $*(y(s))^(s,0,a;)  rfs,  |x| ->  oo, 
X  V  TT  7-00  OS 


(9) 


where  the  integration  is  taken  in  the  clockwise  sense  around  the  profile  sj.  =  0  of  the  airfoil 
S. 

Since  the  noise  sources  are  confined  to  the  neighborhood  of  the  trailing  edge,  we  may 
integrate  by  parts  in  the  final  line  of  (9)  to  express  the  radiation  as  a  Stieltjes  integral 

p(x,u;)  J^^sin2  i  '^{s,0,uj)  d^*,  Ixl  oo,  (10) 

|x|  V  TT  7s 

where  the  integration  is  in  the  clockwise  direction  about  S.  In  §3  it  will  be  convenient  to 
evaluate  the  integral  by  setting  z  =  xi  +  ix2,  and  mapping  the  region  of  the  2:-plane  outside 
the  airfoil  occupied  by  the  fluid  onto  the  upper  half  of  the  complex  ^-plane,  such  that  the 
airfoil  profile  S  transforms  into  the  real  ^-axis  with  the  correspondence  of  boundary  points 
indicated  in  Figure  2b.  But,  the  potential  is  given  in  terms  of  C  by 


$*  =  $*(^)  =  -pReC,  /i>0,  (11) 

where  p  is  a  constant  whose  value  depends  on  the  shape  of  the  trailing  edge  profile.  Thus 
the  acoustic  pressure  may  be  written 

p{x,u)  «  -^^^W^^sin^  'i/'sin(0/2)e*''°l'‘l  f  4^(s(C),0,a;)  d(,  |x|  oo,  (12) 

|X|  V  J~oo 

where  the  path  of  integration  runs  just  above  any  singularities  on  the  real  axis. 
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3.  MODELING  HIGH  FREQUENCY  EDGE  NOISE 


3.1  Conformal  mapping 

The  fluid  region  bounded  internally  by  a  trailing  edge  profile  of  the  type  shown  in  Figure 
2a  is  mapped  into  the  upper  ^-plane  by  means  of  the  transformation 

|  =  7j^[/(C,l)+a/(C.^)],  a>0,  /?>+!,  (13) 

where  a  and  /?  are  numerical  constants,  and 


/(<•« 


+  -\n 

TT 


(c  +  ^)  + 


(14) 


Equation  (13)  maps  the  real  (^-axis  (Figure  2b)  onto  an  airfoil  profile  S  whose  upper 
and  lower  surfaces  coincide  with  X2  =  ±|h  respectively  for  Re  ^  <  — /?  and  Re  >  1. 
This  transformation  divides  the  upper,  rounded  section  of  the  trailing  edge  (between  B 
and  A  in  Figure  2a)  in  two  sections:  the  interval  — /?  <  Re  C  <  —  1  is  the  image  of  the 
section  that  starts  at  B  and  terminates  (at  C  =  —  1)  at  the  top  of  a  vertical  ‘end-face’ 
(oTi  =0,  —  <  0)2  <  —  -f  A)  of  thickness  A,  which  corresponds  to  the  interval 

— 1  <  Re  C  <  1.  The  values  of  a  and  (5  are  calculated  by  prescribing  values  for  the 
geometrical  ratios  i/h  and  A/h.  As  in  Chapter  1,  results  presented  in  this  chapter  are  for 


a  =  600,  0  =  86.9370 


(15) 


which  yield  £/h  =  4  and  A/h  =  0.0074.  The  corresponding  airfoil  profile  is  shown  in  Figure 
2a  (because  A  <C  h  the  end-face  cannot  be  distinguished). 

The  coefficient  n  in  the  definition  (11)  of  is  determined  from  condition  (4),  which  is 
equivalent  to  —iJ,(  ->  —iy/z  as  \z\/h  — >  oo,  and  supplies 


H  = 


h 

7r(l  -f  a) 


(1 


4q; 

+w 


(16) 


3.2  The  frozen  approximation 

In  the  frozen  approximation  of  thin  airfoil  theory  [5  -  7]  the  boundary  layer  turbulence 
is  assumed  to  convect  parallel  to  the  undisturbed  mean  stream  at  a  convection  velocity  Uc, 
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which  depends  weakly  on  frequency  ui.  When  this  is  applied  to  the  airfoil  of  Figure  2a, 
boundary  layer  eddies  translate  as  frozen  distributions  of  vorticity  in  the  Xi-direction  within 
a  boundary  layer  constrained  to  lie  above  the  plate  in  the  X2  >  \h,  so  that  the  separation 
streamline  is  parallel  to  the  undisturbed  mean  flow  velocity,  and  extends  downstream 
parallel  to  the  Xi  axis.  ^  satisfies  Laplace’e  equation  in  the  region  X2  <  5/1,  where  we  can 
write 

/OO  1 

X2  < -h.  (17) 

-OO  2 

It  was  shown  in  Chapter  1  (§3.2)  that  the  pressure  fluctuations  Pi(x,t)  produced  by  the 
boundary  layer  turbulence  in  X2  <  \h  when  the  presence  of  the  airfoil  is  ignored  satisfies 
the  linear  relation  dvi/dt  =  (— l/po)Vpi,  and  it  therefore  follows  by  applying  this  relation 
on  X2  =  5/1  that 

A(k,uj)  =  ^^^^Ps{k,0,u}),  (18) 

2poU) 

Ps{k,k3,w)  being  the  Fourier  transform  of  the  blocked  wall  pressure,  defined  by 

Ps{k,k3,uj)  = -^  f  dxidx^dt,  (19) 

[ZTT)^  J-oc 

where  Ps{xi,xz,  t)  =  2pi{xi,h/2,  xz,  t)  is  the  pressure  that  would  be  exerted  by  the  boundary 
layer  on  a  rigid  wall  at  X2  =  ^h,  i.e.,  on  the  upstream,  flat  section  of  the  upper  surface  of 
the  airfoil.  In  the  frozen  approximation  it  is  assumed  that  measurements  of  Ps{xi,xz,t) 
several  characteristic  hydrodynamic  length  scales  upstream  of  the  trailing  edge  are  sufficient 
to  determine  Ps{k,  kz,io). 


Then,  in  equation  (12)  we  find 

/OO  ^ 

«(s(0,o,i.,)rfc 

-OO 

X(k) 


1 

- -  /  X{k)ps{k,Q,u))  dk, 

2PqLO  J — OO 

C.C.  l^Q-kh/2  J^^Q-ikz 

< 

— c.c.{J(— ^)}  for  k  <  0, 


(20) 

(21) 


where  ‘c.c.’  denotes  the  complex  conjugate.  Thus,  equations  (12)  and  (1)  give  the  far  field 
sound  in  the  form 


p(x,  t) 


fjL  sin2  ipsm{6/2) 
\/27rCo|x| 


1 dkdw,  |x|  00.  (22) 


The  acoustic  pressure  frequency  spectrum  $(x,ci;),  which  satisfies 

(p2(x,t))  =/o°°$(x,a;)da;. 
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the  angle  brackets  (  )  denoting  an  ensemble  average,  is  calculated  from  (22)  by  first 
assuming  that  only  a  finite  spanwise  section  —\L<Xz<\Loi  the  trailing  edge  is  wetted 
by  the  turbulent  flow,  where  L  is  much  larger  than  the  boundary  layer  thickness  5.  For 
statistically  stationary  turbulence 

L 


(ps(fc,0,c<;)p*(fc',0,a;'))  «  —6{ijJ  -  u')5{k  -  A:')P(A:,  0,  w),  L  »  (J, 

/TT 


(23) 


where  the  asterisk  represents  the  complex  conjugate,  and  F{k,  ks^u)  =  F{—k,  k^,  —uj)  is  the 
wavenumber-frequency  spectrum  of  the  blocked  wall  pressure  [8].  Then  equation  (22)  yields 

fj?  sm 'll)  {6/2)  ujL  f°° 


$(x,a;)  = 


27r2|x| 


/OO 

\I{k)\‘^F{k,0,uj)dk,  |xHoo. 

-CX) 


(24) 


To  illustrate  predictions  of  this  formula  the  wall  pressure  spectrum  will  be  approximated 
by  the  Corcos  formula,  which  is  applicable  for  lvS/U  >  1  in  the  immediate  neighborhood  of 
the  ‘convective  ridge’  (where  P{k,k3,u)  attains  a  large  maximum)  [9]: 

A  ^3 


P{k,k3,uj)  =  ^pp{u) 


7r[l  +  il{k  —  uj/Uc)'^]  7r[l  -f 


«  9Uc/uJ,  4  «  lAUc/oj. 


(25) 


The  lengths  and  4  are  frequency  dependent  turbulence  correlation  scales  in  the  Xi-  and 
a:3-directions.  $pp  is  the  point  pressure  frequency  spectrum  which  will  be  approximated  by 
the  following  empirical  formula  (based  on  data  collated  by  Chase  [8]): 

iU/S,)%p{oj)  _  (uS./U)^ 


ap  —  0.12, 


(26) 


where  <5*  ~  5/8  is  the  boundary  layer  displacement  thickness,  and  is  the  friction  velocity. 

The  principal  contribution  to  the  integral  in  (24)  is  from  the  immediate  vicinity  of  the 
spectral  maximum  in  the  convective  ridge,  where  k  ~  w/f/c-  We  therefore  replace  X{k)  by 
I{lo/Uc)  to  obtain 

p? sin'll) six? {6/2)  (jjL  ^  /-oo 


$(x,a;) 


27r2|x| 


^(1)1 


Q.lLMcP?  sin  ij)  sxx?{6/2) 


i{^)\  *„(-). 

where  the  Me  =  Uc/Co  is  the  convection  Mach  number. 


TT^IxP 


oo. 


(27) 
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The  integral  in  (21)  defining  X{k)  must  be  evaluated  numerically,  with  z  =  2:(C) 
determined  by  equation  (13).  Then  the  acoustic  frequency  spectrum  can  be  expressed  in 
the  nondimensional  form 

(t//(i*)$(x,a;)  ^  L  /  o;  \  P  ao{u;5./U) 


[povD 


2\2 


/i^M  (  )  sin'^sin^ 


©(^)K 


U) 

Uc 


[al  +  {u;6jU)^]l 


gq  —  0.035, 
(28) 


where  M  =  U/co  is  the  Mach  number  of  the  main  stream,  and  we  have  introduced  the 
approximation  Uc  ~  0.7U;  //  is  given  by  (16). 


The  solid  curve  in  Figure  3  is  a  nondimensional  plot  of  the  acoustic  pressure  spectrum  for 
the  particular  case  h  =  6,  when  the  thickness  of  the  turbulent  boundary  layer  over  the  upper 
surface  is  just  equal  to  the  airfoil  thickness  h.  At  very  low  frequencies  (when  uh/Uc  1) 
the  radiation  must  be  similar  to  that  from  the  edge  of  a  thin  rigid  half-plane.  In  this  limit 


(I) 


ttUc 


u)h 


fi  \  oj  +  i0^ 

and  the  spectrum  becomes  identical  with  that  for  the  half  plane 

(U/i.)4.(x,u.)  ^  „  (L5.\ ,,, ,  fe\  aM.fU) 


which  is  also  plotted  in  Figure  3 


M 


j  sin  tp  sin^  [  ^  ]  3 

|x|V  VJ  [al  +  {u;6^/Uy]l 


(29) 


(30) 


The  large  difference  at  high  frequencies  between  the  half  plane  and  rounded  edge  spectra 
occurs  because,  in  the  frozen  approximation  the  perpendicular  stand-off  distance  of  the 
separation  streamline  from  the  edge  A  (where  the  high  frequency  noise  is  generated)  is 
equal  to  the  plate  thickness  h,  so  that  the  upwash  velocity  at  the  edge  is  reduced  by  a  factor 
of  order  <C  1  relative  to  the  half  plane. 


3.3  Fully  attached  edge  flow 

When  the  turbulent  boundary  layer  remains  attached  to  the  curved  upper  surface  of  the 
airfoil  right  up  to  the  edge  A,  as  indicated  in  Figure  4,  we  write  the  representation  (17)  for 
4'(a:i,a:2,  w)  on  the  upper  surface  in  terms  of  (s,  Sx)  in  the  form 

^(s,sx,a;)=  /  dk,  s±<0,  C  <  1-  (31) 

A{k,u))  is  determined  by  assuming  that  the  wall  pressure  Fourier  transform  Ps{k,0,oj)  is 
defined  locally  on  the  curved  surface,  with  local  statistical  properties  defined  as  in  (23)  -  (26) 
that  may  vary  over  length  scales  much  larger  than  the  correlation  scale  of  the  turbulence. 


45 


Report  No.  AM  99-003 


Boston  University,  College  of  Engineering 


This  approximation  can  be  used  to  evaluate  the  portion  4'(s(C),  0,a;)  d(^  of  the 
integral  in  the  representation  (12)  of  the  sound  pressure.  On  the  lower  surface  of  the  airfoil 
(i.e.,  for  1  <  C  <  +oo),  we  apply  a  local  plane  wave  approximation  to  (31)  by  writing 

/OO 

sj.  =  0,  C  >  1,  (32) 

-OO 

where  for  a  point  on  the  lower  surface  specified  by  C  (1  <  C  <  +oo)  s'(C)  is  the  value  of 
s  at  the  corresponding  point  on  the  upper  surface  (i.e.,  where  the  perpendicular  from  the 
lower  surface  point  meets  the  upper  surface),  h'{C)  is  the  airfoil  thickness  at  that  point,  and 
d  =  i?(^)  is  the  angle  between  the  sj.-axis  and  the  a:2-direction  at  the  upper  surface  point. 
Thus,  when  A{k,u)  is  approximated  in  terms  of  the  blocked  pressure  as  in  (18),  the  far 
field  sound  is  found  to  be  given  by 


p(x,  t) 


-/isins  ipsm{9/2) 
|x|v'27rCo 


ff  dkdio  f  sgn{k)y/iiops{k,0,Lo)e^'‘^^^^ 

J  J —  OO  j —OO 


X 

where,  for  — co  <  C  <  Ij  C(C)  is  determined  by  the  point  of  intersection  with  the  lower 
surface  of  the  perpendicular  onto  the  latter  from  the  point  C  on  the  upper  surface,  so  that 
1  <  C  <  -f-oo,  and  d  =  ^^(O- 


1  I  /  -ifcV(C)sintf— |fc|/i'(^)costf 


g-ta.(t-|x|/c„)  (33) 


Thus,  forming  the  acoustic  pressure  frequency  spectrum,  as  in  §3.2,  we  obtain  the 
following  modified  form  of  equation  (24) 

//^sin^sin^(^/2)  wL  r®® 


$(x,a;) 


X 


1  + 


27r2|x| 

A. 


—  f  dk  [f  P(A:,0,a;)e" 

Cq  J  —  OO  j  j —  OO 


ifc{s(C)-s(C')} 


ikh'{()  sin 'd—\k\h'{Q  cos d  1  |  |  — ikh'{^')  sin 'd'—\k\h'{(')  cos'd' 

J  r  Uc7 


X  OO. 


dCdC, 

(34) 


We  can  take  account  of  slow  changes  (on  a  scale  of  the  boundary  layer  thickness  5)  in  the 
mean  properties  of  the  boundary  layer  turbulence  with  curvilinear  distance  along  the  airfoil 
by  permitting  P(A:,0,a;)  in  (34)  to  be  a  slowly  varying  function  of  s  =  $(C)  and  s'  =  s(C0- 
The  simplest  way  to  do  this  is  to  use  a  modified  form  of  the  Corcos  approximation  (25) 
in  which  the  point  pressure  frequency  spectrum  ^ppiw)  is  replaced  by  the  geometric  mean 
y/$pp(s,  u)^pp(s',  oj)  of  the  corresponding  spectra  at  s  and  s',  viz. 

P(k,  k3,U>)  =  ^^pp(s,  U)^pp(s',  ^[1  +  £p2]  •  (35) 
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This  would  be  appropriate  if  variations  in  the  convection  velocity  Uc  ^  0.7U  are  small 
compared  with  changes  in  the  structure  of  the  point  pressure  spectrum  as  the  trailing 
edge  is  approached.  It  is  necessary  to  take  the  geometric  mean  (as  opposed  to,  say,  the 
arithmetic  mean)  to  ensure  that  the  right  hand  side  of  (34)  remains  positive  definite. 

Thus,  when  the  wavenumber  integral  in  (34)  is  evaluated  as  before,  by  taking  the  first 
term  in  the  expansion  about  the  convective  peak  at  =  uj/Uc,  we  find 


$(x,a;) 


0.7/x^LMsin^sin^(^/2) 

7r^|x  ^ 


x|  — )■  oo, 


—  f  ^$pp(s,  Lo) 
J—oo 


2  I  ^  I  gjfc/i'(C)sint?-|fc|/i'(C) cost9 


Jks{0 


dC 


(36) 


In  the  numerical  results  to  be  given  below  we  shall  ignore  the  dependence  of  ^pp{s,u>)  on 
s.  In  that  case  the  acoustic  pressure  spectrum  has  the  form  given  previously  in  equation 
(28),  with 


I{k)  =  sgn{k)  C  1  +  (  — ) 

J—oo  \  d(^  I 


Jks{0 


dC 


(37) 


This  integral  must  be  evaluated  numerically.  In  doing  this  it  has  been  assumed  that  the 
separation  streamline  at  the  edge  makes  an  angle  —  |  with  the  positive  direction  of  the 
xi-axis,  which  is  the  same  as  for  potential  flow  past  the  airfoil  when  the  Kutta  condition  is 
applied  at  A  (see  Chapter  1).  In  order  to  achieve  this  smooth  behavior,  we  have  set  ^  | 

at  those  points  very  close  to  A  where  the  geometrical  value  of  the  angle  actually  exceeds  |. 


The  acoustic  pressure  frequency  spectrum  is  given  by  equation  (28).  It  is  plotted  in 
Figure  4  for  /i/5  =  1  and  5.  The  intensity  of  the  sound  decreases  as  h/S  increases,  because 
the  exponential  decay  of  the  second  term  in  the  square  brackets  of  (37)  rapidly  reduces  the 
contribution  from  the  lower  surface.  However,  the  full  effect  of  scattering  by  the  upper 
surface  is  always  present,  so  that  the  overall  reductions  with  increasing  thickness  are  very 
much  less  than  those  predicted  by  the  frozen  approximation. 


3.4  Separated  edge  flow 

When  the  boundary  layer  separates  at  some  intermediate  point  on  the  curved  section  of 
the  trailing  edge  (Figure  5)  the  radiation  can  be  estimated  by  combining  the  approximations 
of  §§3.2,  3.3.  Let  separation  occur  at  xi  =  Xo  {C  =  Co),  the  point  where  the  airfoil 
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thickness  is  ho,  and  assume  that  immediately  after  separation  the  turbulence  convects  in  a 
frozen  pattern  parallel  to  the  undisturbed  mean  stream.  In  the  simplest  approximation  in 
which  the  dependence  of  $pp(s,a;)  on  s  is  again  ignored,  the  acoustic  pressure  spectrum  can 
be  set  in  the  form  (28),  where 


<ic|,  (38) 

where  d{C)  =  ho  —  \h  —  X2{C)  (C  <  1)  is  the  vertical  standoff  distance  of  the  separation 
streamline  from  the  upper  surface  of  the  airfoil,  and  So  =  s(Co)-  The  separation  coordinate 
Xo  and  the  curvilinear  distance  of  the  separation  point  from  the  edge  are  displayed  in  Table  1. 


I{k)  =  sgn(A:)< 


1-f 


A. 


^ikh'iC)  sin  cos  d 


ifcs(C) 


dC 


+ 


f 

JCa 


.-imo 


-1- 


A, 


\k\ho 


^ik(so+xi-Xo) 


ho/h 

Xo/h 

curvilinear  distance/h 

1 

-4 

4.19 

0.5 

-0.97 

1.11 

0.2 

-0.24 

0.32 

0.1 

-0.08 

0.13 

Table  1.  Separation  coordinates  for  a  —  600,  P  =  86.9370. 

Typical  acoustic  pressure  frequency  spectra,  calculated  from  the  general  formula  (28), 
are  illustrated  in  Figure  5  for  the  cases  of  Table  1  when  S  =  h.  The  contributions  at  high 
frequencies  decay  very  rapidly  as  the  separation  point  moves  upstream  from  the  edge,  for 
the  reasons  discussed  in  §3.2.  The  case  ho/h  =  1  coincides  with  the  frozen  approximation 
of  §3.2. 

These  results  (and  the  frozen  approximation  of  §3.2)  take  no  account  of  the  sound 
generated  by  turbulence  in  the  ‘dead  water’  region  bounded  by  the  separation  streamline 
and  the  curved  section  of  the  upper  surface  between  the  separation  point  and  the  edge. 
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4.  CONCLUSION 

The  intensity  of  the  high  frequency  sound  generated  by  turbulent  flow  past  the  trailing 
edge  of  an  airfoil  is  significantly  reduced  if  separation  occurs  upstream  of  the  geometrical 
edge.  Separated  flow  over  a  curved  trailing  edge  profile  (Figure  5)  can  be  modeled  by 
assuming  that  after  separation  boundary  layer  eddies  are  convected  along  a  path  parallel 
to  the  undisturbed  mean  stream,  so  that  their  distance  of  closest  approach  to  the  sharp 
trailing  edge  is  just  equal  to  the  airfoil  thickness  at  the  separation  point.  The  strength 
of  the  unsteady  interaction  responsible  for  sound  generation  then  decreases  exponentially 
with  increasing  frequency,  and  predicted  levels  of  edge  noise  are  significantly  smaller  than 
estimates  made  by  modelling  the  interaction  in  terms  of  a  rigid  half  plane.  This  conclusion 
ignores  possible  contributions  to  sound  generation  from  turbulence  in  the  recirculating  zone 
between  the  separation  streamline  and  the  curved  section  of  the  upper  surface  between  the 
separation  point  and  the  edge,  which  are  probably  significant  only  at  much  lower  frequencies 
[!]• 

When  the  mean  flow  remains  attached  right  up  to  the  trailing  edge,  boundary  layer 
turbulence  on  the  suction  side  of  the  airfoil  always  interacts  strongly  with  the  edge. 
However,  because  of  the  finite  thickness  of  the  airfoil,  the  influence  of  the  (turbulence  free) 
lower  surface  on  the  radiated  intensity  continues  to  decrease  exponentially  fast  at  high 
frequencies.  In  consequence,  most  of  the  high  frequency  radiation  is  associated  with  the 
flow  interacting  with  the  upper  surface.  The  sound  levels  are  smaller  than  for  the  half  plane 
model,  but  the  differences  at  higher  frequencies  are  more  modest,  typically  being  of  the 
order  5  to  10  dB,  the  precise  values  being  dependent  on  the  ratio  of  the  boundary  layer 
thickness  to  the  mean  airfoil  thickness. 
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Figure  1.  Schematic  low  Mach  number  turbulent  boundary  layer  flow  over  a  rounded 
trailing  edge. 
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Figure  2.  (a)  Configuration  of  the  trailing  edge  and  the  coordinate  system. 

(b)  Mapping  the  region  outside  the  airfoil  (in  the  plane  of  ^  =  Xi+ix2) 
onto  the  upper  half  of  the  C-plane. 
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Figure  3.  Frozen  approximation  to  the  edge  noise  spectrum 

(i7/(5*)$(x,a;)/[ao(/Oo'y*)^M(L(5*/|xp)  sin?/)sin^(0/2)] 

for  the  rounded  edge  of  Figure  2a  of  thickness  h  =  5  when  I  —  Ah, 

and  for  the  rigid  half  plane  (/i  =  0). 
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SEPARATED  EDGE  FLOW 


Figure  5.  Predicted  edge  noise  spectrum 

(f//J*)$(x,w)/[oo(pot;2)2M(L5*/|xp)sint/)sin^(0/2)] 

for  the  rounded  edge  of  Figure  2a  when  h/(5  =  1,  £  =  4/i,  and  for  the 

separated  edge  flows  of  Table  1. 
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